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I. INTRODUCTION 

In the infinite quark mass limit, heavy quark spin-flavor symmetry (HQSFS) implies that the 
dynamics involving heavy quarks are independent of their spin or flavor. In this way, charm and 
bottom spectra can be related, up to corrections suppressed by 1 /mg with ttiq the heavy quark 
mass. It should also be possible, in a given heavy flavor sector, to relate states with different spins. 
These relations are very useful in the study of composites that mix heavy and light quarks. In this 
work, we focus on hadronic molecular states composed by a heavy-light meson and a heavy-light 
antinreson (p(*)p(*), P = D,B). These molecular states were predicted in the mid 70s [1, 2]. So 
far, the best experimental candidate to fit this molecular description is the X (3872) resonance, first 
observed by the Belle Collaboration in 2003 [3], that can be thought as a DD* bound state with 
J PC = 1 ++ (quantum numbers confirmed later on in Ref. [4]). Since then, many other new XYZ 
states which are good candidates to be exotic hadrons have been experimentally observed [5, 6]. 

Within the molecular description of the X(3872), the existence of a X 2 [J pc = 2 ++ ] s-wave 
D*D* bound state was predicted in the effective field theory (EFT) approach of Refs. [I, 8]. As 
a result of the heavy quark spin symmetry (HQSS), the binding energy of the X 2 resonance was 
found to be similar to that of the A(3872), i.e., 

Mx 2 ~ A/x( 3872 ) ~ Md* — Md ~ 140 MeV. (1) 

The existence of such a state was also suggested in Refs. [9-14]. Both the X(3872) and the X 2 , to 
be denoted by ^2(4013) in what follows, have partners in the bottom sector [15], 1 which we will 
call Xb and X^, respectively, with masses approximately related by 

Mx b2 — XIx b ~ Mb* — Mb ~ 46 MeV. (2) 

It is worthwhile to mention that states with 2 ++ quantum numbers exist as well as spin partners 
of the 1 ++ states in the spectra of the conventional heavy quarkonia and tetraquarks. However, 
the mass splittings would only accidentally be the same as the fine splitting between the vector 
and pseudoscalar charmed mesons, see Eq. (I). 2 For instance, the mass splitting between the first 
radially excited charmonia with 2 ++ and 1 ++ in the well-known Godfrey-Isgur quark model is 
30 MeV [17], which is much smaller than the value in Eq. (1). In a quark model calculation with 
screened potential, the 2 ++ — 1 ++ mass splitting for the 2 P charmonia is around 40 MeV [18]. As 

1 In Ref. [15], the bottom and charm sectors are connected by assuming the bare couplings in the interaction La- 
grangian, see Eq. (A6), to be independent of the heavy quark mass. This assumption will also be used throughout 
this work. 

2 Were these states due to threshold cusps, the splittings would be the same as those of the hadronic molecules. 
However, it was shown in Ref. [16] that narrow near threshold peaks in the elastic channel cannot be produced by 
threshold cusps. 
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for the tetraquark states, the corresponding mass splitting predicted in Ref. [19] is 70 MeV, which 
is again much smaller than Mp* — Mp. Notice that it is generally believed that the x C 2 ( 2 P) has 
been discovered [20, 21], and its mass is much lower than 2 Mp*. Therefore, we conclude that a 
possible discovery of a 2 ++ charmonium-like state with a mass around 4013 MeV as a consequence 
of HQSS [ 22 ] would provide a strong support for the interpretation that the V(3872) is dominantly 
a DD* hadronic molecule. It is thus very important to search for such a tensor resonance, as well 
as the bottom analogues, in various experiments and in lattice QCD (LQCD) simulations. 

Some exotic hidden charm sectors on the lattice have been recently studied [23-27], and evidence 
for the A"(3872) from DD* scattering on the lattice has been found [24], while the quark mass 
dependence of the X(3872) binding energy was discussed in Refs. [28, 29]. The 2 ++ sector has 
not been exhaustively addressed yet in LQCD, though a state with these quantum numbers and 
a mass of (m Vc + 1041 ± 12) MeV= (4025 ± 12) MeV, close to the value predicted in Refs. [7, 8 ], 
was reported in Ref. [23]. The simulation used dynamical fermions, novel computational techniques 
and the variational method with a large basis of operators. The calculations were performed on 
two lattice volumes with pion mass ~ 400 MeV. There exists also a feasibility study [30] of future 
LQCD simulations, where the EFT approach of Refs. [7, 8 ] was formulated in a finite box. 

On the other hand, despite the theoretical predictions on their existence, none of these hypo¬ 
thetical particles has been observed so far. Nevertheless, they are being and will be searched for 
in current and future experiments such as BESIII, LHCb, CMS, Belle-II and PANDA. It is thus of 
paramount importance to provide theoretical estimates on their production rates in various experi¬ 
ments, as well as the dominant decay modes and widths . 3 The production of these states in hadron 
colliders and electron-positron collisions has been studied in Refs. [31, 32]. In this work, we will 
investigate the dominant decay modes of the spin-2 partners of the A(3872), i.e. the A*2(4013) and 
Xb 2 , and provide an estimate of their decay widths. 

Besides, we will also discuss the radiative X 2 —» DD *7 and X 1,2 —» BB *7 transitions. These 
decay modes are more sensitive to the long-distance structure of the resonances and might provide 
valuable details on their wave-functions. The situation is similar to that of the V(3872) —>■ D°D°tt° 
decay studied in Refs. [33, 34]. Also here, the widths will be affected by the DD* or BB* final 
state interaction (FSI). FSI effects are expected to be large because they should be enhanced 
by the presence of the isovector Z c (3900) and 2^(10610) resonances located near the DD* and 
BB* thresholds, respectively. Besides, FSI corrections will be also sensitive to the negative C- 
parity isoscalar DD* or BB* interaction. Eventually, precise measurements of these radiative 

3 If a resonance is too broad, say F > 200 MeV, it would be very difficult to be identified since it is highly nontrivial 
to distinguish the signal for a broad resonance from various backgrounds. 
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decay widths might provide valuable information on the interaction strength in this sector, which 
would be important in understanding the p(*)p(*) system and other exotic systems related to it 
through heavy quark symmetries [15, 35]. 

The structure of the paper is as follows. First in Sect. II, we briefly discuss the relation of 
the charm and bottom 2 ++ states with the X(3872) resonance, and in Sect. Ill we present our 
predictions for the X 2 —>• DD , DD* hadron decays and the X^ 2 BB,BB* ones in the bottom 
sector. In Sect. IV, the X 2 and X^ radiative decays are investigated, paying special attention to the 
loop mechanisms responsible for the FSI contributions. The conclusions of this work are outlined 
in Sect. V and in addition, there are three Appendices. In the first one (Appendix A), we collect 
different heavy meson Lagrangians used through this work, while the validity of the perturbative 
treatment of the DD for the X 2 is discussed in the second one (Appendix B). Finally, in Appendix C, 
we give some details on the evaluation of different three-point loop functions that appear in the 
computation of the hadronic and radiative decays. 


II. HQSFS, THE X(3872) RESONANCE AND THE CHARM AND BOTTOM X 2 STATES 


A. X (3872) 


As mentioned, we start assuming the X (3872) to be a positive C-parity DD* bound state, with 
quantum numbers J pc = 1 ++ . At very low energies, the leading order (LO) interaction between 
pseudoscalar and vector charmed (D°, D + , D*°, D * + ) and anti-charmed (Z)°, D ~, D*°, D *~) mesons 
can be described just in terms of a contact-range potential, which is constrained by HQSS [7, 8, 15]. 
Pion exchange and particle coupled-channel 4 effects turn out to be sub-leading [7, 36]. For the case 
of the X(3872), isospin breaking is important [37] as this bound state is especially shallow. The 
energy gap between the D°D*° and D + D*~ channels is around 8 MeV, which is much larger than 
the A(3872) binding energy with respect to the D°D*° threshold. As a consequence, the neutral 
(D°D*°) and charged ( D + D *~) channels should be treated independently. The coupled-channel 5 
contact potential in the 1 ++ sector is given by [8] (see also Appendix A 1) 

1 ( Cox + C\x Cox — C\x 


V : 


X(3872) 


Cox — C\x Cox + Cix 


(3) 


with Cox and C\x low energy constants (LECs) that need to be fixed from some input. This 
interaction is used as kernel of the Lippmann-Schwinger equation (LSE) in the coupled channel 

4 We do not refer to charge channels, but rather to the mixing among the DD, DD *, D*D* pairs in a given IJC 
(isospin, spin and charge conjugation) sector. 

5 Actually, positive C-parity combinations in both the neutral D° D*° and charged D + D*~ channels are being 


considered. 
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space in the l +_l “ sector, 

T(E-p',p) = V(p',p) + 


d 3 q 


V(p',q)- 


— T(E;q,p) 


(4) 


(2vr) 3 ^ E- q 2 /2p l2 -M l -M 2 + ie 

with M± and M 2 the masses of the involved mesons, fiy 2 = Af] -1 + M^ 3 , E the center of mass (c.m.) 
energy of the system and p ( p') the initial (hnal) relative three momentum of the DD* pair in the 
c.m. frame. The used normalization is such that above threshold [E > (Mi + M 2 )], the single 


channel elastic unitary condition is ImT _1 (.E) = p\ 2 k/{2'K) with k = ^ J2[i\ 2 (E — M x — M 2 ). The 
discussion is similar for any other J PC sector. Due to the use of contact interactions, the LSE 
shows an ill-defined ultraviolet (UV) behavior, and requires a regularization and renormalization 
procedure. We employ a standard Gaussian regulator (see, e.g. [38]) 

<: p\V\p ') = C IX e^ 2/A2 e~v' 2 / K \ (5) 


with Cjx any of the LECs of Eq. (3) in the case of the X(3872), or the relevant ones for any other 
J PC sector. We take cutoff values A = 0.5 — 1 GeV [7, 8], where the range is chosen such that A 
will be bigger than the wave number of the states, but at the same time it will be small enough 
to preserve HQSS and prevent that the theory might become sensitive to the specific details of 
short-distance dynamics. 6 The dependence of the results on the cutoff, when it varies within this 
window, provides a rough estimate of the expected size of sub-leading corrections. Bound states 
correspond to poles of the T-matrix below threshold on the real axis in the first Riemann sheet (RS) 
of the complex energy, while the residues at the pole give the s-wave couplings of the state to each 
channel ( D°D*° and D + D*~ in the case of the X(3872) resonance 7 ). 

The LECs Cox and C\ X can in principle be determined [8] from M X (;is 72 ) = (3871.69 ± 0.17) 
MeV (mass average quoted by the PDG [5]) and the isospin violating ratio of the decay amplitudes 
for the X(3872) — > J/i/jtttt and X(3872) — > J/ipTnnr, R X ( 3872 ) = 0-26 ± 0.07 [39]. The ratio is given 
by (see Eq. (80) of Ref. [8]) 


R-X(3872) 


dh, + 


( 6 ) 


where ^ n>c give the average of the neutral and charged wave function components in the vicinity of 

the origin, and are related to the LECs introduced in Eq. (3) by [37] 

'kn _ 1 -G 2 (Cox + C xx ) /2 _ G x (C ox - C 1X ) /2 

T c C 2 (Cox - C\x) /2 1 - Ci (Cox + C 1X )/2 U 


6 However, as will be shown later on, the situation is more complicated in the two-body d-wave hadronic decays. 

7 For instance, in the case of the X(3872), we have 

i g °y = ,e ■ Mx(3872>i x ru(E) ■ 

goQc = lim [E - M x (3872)] x T 12 (E) , 

E—tMx ( 3872 ) 

where Tij are the matrix elements of the T-matrix solution of the UV regularized LSE. 
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where the neutral and charged loop functions, G i = G d oq*o = G D *ojjo,G 2 = G D + D *~ = G D *+ D ~, 
are defined in Eq. (A12), and should be evaluated at the X(3872) pole mass. We use m D o = 
(1864.84 ± 0.07) MeV, m D + = (1869.61 ± 0.10) MeV, m D *o = (2006.96 ± 0.10) MeV and m D *+ = 
(2010.26 ± 0.07) MeV [5]. Note that mjy o + o = (3871.80 ± 0.12) MeV, and the uncertainty 
in the value of this lowest threshold affects the precision of the X(3872) binding energy. We have 
taken into account this effect by adding in quadratures the PDG error of the X(3872) mass and 
that of the neutral channel threshold and assign this new error to the mass of the resonance, that 
now reads M X (ss7 2 ) = (3871.69 ± 0.21) MeV. For the LECs, we obtain: 

Cox = -1.70l°;°? (-0.731™) fm 2 , C iX = -0.09lg;^ (-0.38t°;} 2 ) fm 2 , (8) 

for A = 0.5(1.0) GeV. Errors, at the 68% confidence level (CL), have been obtained from a 
Monte Carlo (MC) simulation assuming uncorrelated Gaussian distributions for the two inputs 
(Mw 38 72 ), T?x( 3872 ))■ I n the simulation, we have rejected MC samples for which the A(3872) 
turned out to be unbound, since the scheme of Ref. [8] only allows to determine the properties 
of the resonance when it is bound. 


B. X2(4013): J pc = 2++, charm sector 


HQSS predicts that the s-wave D*D* interaction in the 2 ++ sector is, up to corrections sup¬ 
pressed by the charm quark mass, identical to that in the X(3872) sector (1 ++ ) and given by Eq. (3) 
[7, 8]. Thus, in the 2 ++ sector, the potential in the ( D*°D *°), ( D* + D*~) coupled channel space 
reads (see also Appendix A 1) 


V 2 ++ = 


+ 0(q/m c ), 


(9) 


1 I Cox + Cix Cox - C\x 
Cox — Cix Cox + C±x 

with the same structure and involving the same LECs that in the X(3872) channel. Besides, in 
the above equation m c ~ 1.5 GeV is the charm quark mass and q ~ Aqcd> a scale related to the 
light degrees of freedom. Taking Aqcd ~ 300 MeV [5], corrections of the order of 20% to the 
interaction predicted by HQSS cannot be discarded 8 . Nevertheless, it seems natural to expect a 
2 ++ D*D* loosely bound state (X 2 ), the HQSS partner of the X(3872), and located in the vicinity 
of the D*°D*° threshold (~ 4014 MeV) [7, 8, 15]. This is illustrated in the X 2 binding energy 
distributions depicted in Fig. 1. Neglecting the 0(q/m c ) corrections to the LECs, and using those 


8 The two type of uncertainties (HQSS subleading corrections and errors inherited from the inputs) affecting the 
determination of the LECs in the 2 ++ sector are statistically uncorrelated, and should be accordingly added up. 
To that end, we use MC techniques as explained in the caption of Fig. 1. 
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obtained from the A%3872) resonance, we find a clear signal (blue histograms) of a weakly bound 
state. However, the case is less robust when the latter corrections are taken into account. Thus, 
because of the additional 20% HQSS uncertainty, the area below the red shaded A = 0.5 GeV (A = 1 
GeV) histogram is only 0.77 (0.68). This means that approximately a 23% (32%) of A(3872) events 
[{M X ( 3872 ) >( 3872 )) MC samples] do not produce a X 2 pole in the first RS, since the strength of 
the resulting interaction in the 2 ++ sector would not be attractive enough to bind the state, though 
a virtual state in the second RS will be generated instead. Given the existence of the A(3872) as 
a DD* molecule, if the X 2 resonance exists, we would expect its mass (binding energy) to lie most 
likely in the interval [2m£)*o,4006 MeV] ([0,8] MeV), as shown in Fig. 1. Note that the discussion in 
Ref. [15] was simpler, because there we worked in the isospin symmetric limit and used the averaged 
masses of the heavy mesons, which are larger than those of the physical D° and D*° mesons. 



M X ( 4013 ) (MeV) 


M X ( 4013) (MeV) 


FIG. 1. X 2 binding energy histograms obtained from the interaction of Eq. (9) using LECs distributions 
determined from the A(3872) resonance inputs (blue) or using LECs distributions additionally modified to 
account for the HQSS systematic error (red). Left and right plots correspond to UV cutoffs of 0.5 and 
1 GeV, respectively. MC sample (Cox + C±x,Cox — Cix) pairs fitted to the input (M X ( 3872 )>( 3872 )) 
distributions are first generated and are used to evaluate the X 2 mass. A(3872) mass trials above threshold 
are rejected. To evaluate the red shaded histogram, and to account for the HQSS 20% uncertainty in the 
2 ++ interaction, each of the members of any MC sample (Cox + Cix,Cox — Cix) pair is multiplied by 
independent N(/i = 1 ,a = 0.2) Gaussian distributed random quantities r±. 


For later use, we also need the couplings of the X 2 to its neutral (D*°D*°) and charged (D* + D*~) 
components, g and g ^ 2 , respectively. They turn out to be slightly different because the X 2 
resonance is an admixture of isospin 0 and 1, since its binding energy is much smaller than the 
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energy difference between the two thresholds [8, 37]. Considering the HQSS uncertainties, we find: 

10 2 <7(f 2 = 1.4i£i (1.5±Ji) MeV _1//2 , (10) 

10V 2 = l-5±i;| (1.3±5;1) MeV" 1 / 2 , (11) 

for A = 0.5(1.0) GeV. 


C. Xb 2 - J pc — 2++, bottom sector 

Owing to the heavy flavor symmetry, the LO 2 ++ B*B* interaction is given by Eq. (9) as well, 
and thus we should also expect a 2 ++ B*B* bound state ( X fc 2 ), the HQSFS partner of the X(3872), 
located close to the B*B* threshold (~ 10650 MeV) [15]. The X{, 2 binding energy distributions are 
shown in Fig. 2 for the two UV cutoffs employed in this work. We have used the same masses for 
the neutral and charged mesons, m B = (m^o + mg+)/2 = 5279.42 MeV and me* = 5325.2 MeV. 
Note that, according to the PDG [5], — m B o) — ( m B *+ — m B +)\ < 6 MeV CL=95.0%, and 

m B o — m B + = (0.32 ± 0.06) MeV, from where we might expect isospin breaking effects for the B* 
mesons to be significantly smaller than in the charm sector. In Ref. [15] we found that the binding 
energy of the X& 2 state is significantly larger than that of its counterpart in the charm sector (V 2 ), 
around a few tens of MeV. Thus we do not expect any significant isospin breaking effects and the 
Xb 2 resonance would be a pure isoscalar (I = 0) state. 

As can be seen in Fig. 2, in this case we have a robust prediction even when HQSS uncertainties 
(20%) are taken into account. We obtain the mass and the coupling from the residue at the pole 
for A = 0.5 (1.0) GeV: 9 

E Xb2 = 10631%g (105941 22 ) MeV, 10V' 62 = 5.9±^ (6.4+ 2 ;®) MeV~ 1/2 . (12) 

This bound state, being isoscalar, equally couples to the neutral and charged components and, 
therefore: gQ b2 = g^ b2 = g Xb2 . Our predictions in Eq. (12), both for the mass and the B*B* 
coupling of the resonance show some dependence on the UV cutoff, which is to some extent di¬ 
minished when HQSS uncertainties are taken into account. Nevertheless, this A dependence might 
hint to non-negligible sub-leading corrections (among others, pion exchange and coupled channel 
effects [7], which can be larger here than in the charm sector due to the larger binding energy and 
larger meson masses). We will compute the decay widths for both UV regulators, and the spread 
of results will account for this source of uncertainty. 


There appear small differences in the central value of the resonance mass with respect to value quoted in [15] due 
to small differences in the used hadron masses. 
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B Xbl (MeV) 


B Xa (MeV) 


40 35 30 25 20 15 10 5 120 100 80 60 40 20 0 



FIG. 2. Same as in Fig. 1 but in the bottom sector. To better appreciate the distribution details, the A = 
0.5 (1) GeV red histogram, which includes the 20% HQSS error, has been multiplied by a factor of 5 (10). 


III. THE HADRONIC X 2 AND X b2 DECAYS 


The quantum numbers, J pc = 2 ++ , of these resonances constrain their possible decay channels. 
In this work, for hadronic decays we only consider the decays into two heavy hadrons: X 2 —» DD 
and X 2 — > DD*(D*D), and the analogous processes X b2 — > BB and X b2 —> BB*(B*B). We expect 
that these d- wave decay modes should largely saturate the widths of these states. Because the DD 
couples in a d- wave to the 2 ++ system, its contribution to the mass renormalization of the X 2 is 
of higher order (see Appendix B). We thus did not include the DD as a coupled channel in the 
T-matrix, but treat it perturbatively. This means that the transitions are mediated by the exchange 
of a pion. The relevant -n -p(*)p(*) vertices are taken from the LO Lagrangian of heavy meson chiral 
perturbation theory [40-43] (see Appendix A 2). At LO, besides the pion decay constant, = 92.2 
MeV, there appears only one additional D*Dir coupling ( g ). We take g = 0.570 ± 0.006 as inferred 
from the new value of T = (83.4 ± 1.8) keV for the D* + decay width quoted by the PDG [5]. This 
is mostly determined by the recent BABAR Collaboration measurement [44] of this width, which 
is approximately a factor 12 times more precise than the previous value, T = (96 ± 4 ± 22) keV by 
the CLEO Collaboration [45]. Thus, we end up with an uncertainty of the order 1% for g. Though 
the hadronic X 2 and X b2 widths evaluated in this section will be proportional to g 4 , this source of 
error (~ 4%) will be much smaller than others and it will be ignored in what follows. 
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without pion-exchange FF 

with pion-exchange FF 


A = 0.5 GeV 

A = 1 GeV 

A = 0.5 GeV 

A = 1 GeV 

T(X 2 D+D-) [MeV] 

Q 0+3.4 

O.O-i 4 

n q+7.9 
i . 0 _ 2 .! 

U.O—02 

O.8+.2 

T(A 2 D°D°) [MeV] 

9 ' 7 + 3.1 

Z. 1 _! 2 

c 7+7.8 
°- * - 1.8 

n 4+ 05 
u -^-0.2 

0 .6+1 

T(X 2 -i D+D*~) [MeV] 

2 4 + 2 - 1 

4 4+ 31 

0 7 +0 - 6 

u - ‘ -0.3 

1 o +0 - 5 
1 , u - 0.2 

T(A 2 -d- D°D*°) [MeV] 

2 0 +2 1 
z - u -0.9 

0 c+3.5 
°-°- 1.0 

U *°- 0.2 

Q 7 +0.5 
u - * - 0.2 


TABLE I. ^2(4013) —> DD,DD* decay widths using different UV Gaussian regulators for the D*D*X 2 
form factor and with/without including a pion-exchange vertex form factor (FF) in each of the D*Dn and 
D*D*ir vertices in the three-point loop function. The decay widths of the X 2 (4013) —> DD* modes are the 
same thanks to C-parity. Uncertainties are obtained from a Monte Carlo simulation using the X 2 binding 
energy histograms displayed in Fig. 1 (red shaded) and the g* 2 and g* 2 couplings given in Eqs. (10) and 
(11). Note that the procedure takes into account 20% HQSS uncertainties and the correlations between X 2 
masses (binding energies) and g q 2 and g* 2 couplings. Errors on the widths provide 68 % CL intervals. 

A. Charm decays 

1. A 2 (4013) DD 


We will first consider the ^2(4013) —>• D + D~(D°D°) decay, which proceeds through the 
Feynman diagrams depicted in Fig. 3. We treat charm mesons non-relativistically, and neglect 
Pd* d* / m D* terms and the temporal components in the D*,D* propagators. We obtain for the 
^(4013) —>• D + D~ process, in the resonance rest frame and with q and k the 4-momenta of the D 
and D final mesons (q = — k, q° + k° = My 2 ), 


iT(\) D+D 

Ng 2 


f 2 

J 7T 




_ Ng 2 

~ 1 f2 
J 7T 


e»j(A) |.+ 2 J 

/ 


d 4 l 


V V 


(27r) 4 [(l + q) 2 - m? D , + + is] [(k - l) 2 - + ie] ( l 2 - m 2 0 + ie) 

V V 


( 2 + [{l + q) 2 — o + ie] [k — l) 2 — + ie ( l 2 - m 2 _ + ie) 

0) [dc 2 JV { m D *+, 0 ; Mx 2 , (f ) + 2 g* 2 (m D * 0 , m r ; Mx 2 , + )} , 


(13) 


where e t j (A) is the symmetric spin-2 tensor with A denoting the polarization of the X 2 state 10 and 
N = \J 8Mx 2 m^, (y/TriDiTiD *)" accounts for the normalization of the heavy meson fields 11 and some 


10 For instance, it could be given by 

eij(X) = (1,1, 2|Ai, A 2 , A)e\, + 

Ai,A2 

with to = ( 0 , 0 , 1 ), and e±i = T(l, 0 )/\/ 2 , and (j 1 1 j 2 ,j\mi, m 2 , m) a Clebsch-Gordan coefficient. 

11 We use a non-relativistic normalization for the heavy mesons, which differs from the traditional relativistic one by 

a factor V Mh ■ 
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additional factors needed when the couplings g* q, as determined from the residues at the pole of 
the EFT T-matrix, are used for the X 2 D*D* vertex. For the neutral and charged pion masses, we 
have used the values quoted by the PDG [5] and heavy meson isospin averaged masses to compute 
N. Besides, P J is a three-point loop function, the detailed evaluation of which is relegated to the 
Appendix C l. 12 The loop is seemingly logarithmically divergent. However, since the X 2 polarization 
is traceless, the divergent part which comes with a Kronecker delta does not contribute. This is 
because the decay occurs in a d- wave, thus the loop momentum is converted to external momenta, 
and the remaining part of the integral is convergent. Nevertheless, we will include two different 
form factors in the computation of the three-point loop function. One is inherited from the UV 
regularization/renormalization procedure sketched in Eq. (5) and employed to make the LSE T- 
rnatrix finite. In addition, we will include a second form factor to account for the large virtuality 
of the pion in the loop. We will discuss this at length below and in Appendix C 1. 



FIG. 3. Feynman diagrams for the X2(4013) —> D + D decay. Diagrams for the X2(4013) —> D°D° transition 
are similar, with the appropriate changes of the exchanged pion charges. 


Analogously, the A2(4013) —> D°D° amplitude is, 


Ng 2 


-iT{\) D o D o = A) 1 2g* 2 P ] (m D *+, m 7r +; Mx 2 , (f ) + g* 2 1 13 (m rr o , m T o; M X2 , <f 


(14) 


The two-body decay width in the X 2 rest-frame reads [5]: 


dT a 1 


dtt(q) 5 


E ih a ) 


al 327 t 2 M£ 2 ' 


= D + D~, D°D° . 


(15) 


The sum over the X 2 polarizations can be easily done in the c.m. frame, 


e mn(A)ejj(A) 


1 


dmidnj T Onidmj d rnn d t j 


, m,n,i,j = 1,2,3 . 


(16) 


As discussed in Appendix C 1, the three-point loop function has a tensor structure of the type 


I ij (q)=W 2 )q i q j +h(q 2 )S ij \q\ 


(17) 


12 In the computation of P\ we are consistent with the former approximations, and we use non-relativistic charm 
meson propagators. 
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The /] term carries the UV divergence, which however does not contribute to the width, because 
it vanishes after the contraction with the traceless spin-2 polarization tensor, as mentioned above. 
Therefore, only the Jo term is relevant, which is free of UV divergences. Moreover, the contraction 
of 1^(q) I mn (q) with Yl\ e mn(X)e*j(\) [Eq. (16)] leads to a factor of 2<f 4 /3. Thus, the integration 
over the solid angle dfl(q) trivially gives 47r, and the width scales like \q | 5 as expected for a d- wave 
process. 

Our predictions for the V2(4013) —>• D + D~,D°D° decays are compiled in Table I. If we look 
at the first two columns of results in the table, we find widths of the order of a few MeV, with 
asymmetric errors that favour larger values. This is mostly due to the similar asymmetry of the 
uncertainties quoted for the g^ 2 and g* 2 couplings in Eqs. (10) and (11). 

Our scheme is based on a low-energy EFT, in which the momenta should be smaller than a hard 
scale which serves as a momentum cutoff [see Eq. (5)]. The high-momentum modes are out of control 
in the low-energy EFT. Therefore in the computation of the width, we include a Gaussian regulator 
at the D*D*X 2 vertex, as discussed in Eq. (C6). The cutoff should be the same as the one used in 
generating the X 2 as it is related to the same unitary cut in the D*D* system. In Fig. 4, we display, 
as an example, the dependence of the Jo {m D *o , m n o ; Mx 2 , q ^= mp 2 = rn D o ]) integrand [see 
Eqs. (C5)-(C7)] on the pion loop momentum. In the left plot we see that in spite of including the 
Gaussian X 2 form factor, large momenta above 1 GeV provide a sizable contribution to the integral 
(~ 14%, 30% and 45% for A = 0.5 GeV, 1 GeV and 00 , respectively), which is an unwanted feature 
within the low-energy EFT scheme and signals a sizeable short-distance contribution. Indeed, the 
momentum of the exchanged pion, peaks at around 750 MeV, which is a somehow large value in the 
sense that the hard scale for the chiral expansion which controls the pionic coupling is A x ~ 1 GeV. 
We see that below the peak, the curves for both cutoff values are very close to each other, and they 
are also close to the curve corresponding to the case without any regulator. This is the region where 
the low-energy expansion works and thus model-independent conclusions can be made. The curves 
start deviating from each another after the peak, that is in the region with a pion momentum > A x . 
Because the loop integrals are not completely dominated by momentum modes well below A x , the 
widths of interest will bear an appreciable systematic uncertainty. This is reflected in the fact that 
the widths in the second column in Table I are larger than those in the first column by a factor 
around 2. 13 

On the other hand, the fact that the pion could be quite far off-shell should be reflected in the 
D*Dir vertex, which should be corrected, similarly as it is done in the case of the NNtt one. Thus, 

13 Note that the coupling constants obtained with both cutoffs are similar, see Eqs. (10) and (11), and thus the 
difference should come mainly from the loop integration. 



13 


to give an estimate of the hadronic decay widths, in the spirit of the Bonn potential [46], we have 
included a monopole pion-exchange vertex form factor, with a hadron scale of the order of 1 GeV, in 
each of the D*Dn vertices [Eq. (C7)]. Its effect on the internal pion momentum dependence of the 
Io(m D *o, m^o] Mx 2 , = mp 2 = m^o]) integrand is shown in the right plot of Fig. 4. The large 

pion momenta contribution (\l\ > 1 GeV), which is not reliable in a low-energy EFT calculation, is 
reduced now to 6.5%, 13% and 16% for A = 0.5 GeV, 1 GeV and oo, respectively. This makes also 
more appropriate the non-relativistic treatment of the charmed mesons adopted here. Besides, the 
dependence of the width on the UV Gaussian regulator is significantly softer, though the widths 
are further reduced by almost an order of magnitude. 

We believe that the most realistic estimates are those obtained with the inclusion of the pion- 
exchange form factor and the spread of results compiled in Table I give a conservative estimate 
of the systematic uncertainties, beyond the mere existence of the V2(4013) state, as discussed in 
Sect. IIB. We remind here that because of the additional 20% HQSS uncertainty, approximately a 
23% (32%) for A = 0.5 GeV (A = 1 GeV) of the V(3872) events \[Mx{ 3872 ) 5 -Rx( 3872 )) MC samples] 
do not produce the X 2 as a bound state pole, since the strength of the resulting interaction in the 
2 ++ sector is not attractive enough to bind the D*D*. 

Nevertheless, assuming the existence of the X 2 state, and in view of the results given in Table I, 
we estimate the X 2 —>• DD partial width (including both the charged and neutral channels) to be 


T(X 2 -> DD) = (1.2 ± ^03, tH) MeV, 

sys (A) 


(18) 


where the first error accounts for the dependence on the UV Gaussian regulator used in the D*D*X 2 
vertex, while the second one is obtained from the uncertainties given in Table I. This latter error 
includes both some additional systematic (HQSS violations) and statistical (V(3872) input used 
to fix the properties of the X 2 resonance) uncertainties. Notice that, as discussed above, the 
calculation is probably already beyond the valid range of the EFT due to the large contribution 
of high-momentum modes. We thus have adopted a more phenomenological strategy and used the 
pion-exchange form factor with a cutoff of 1 GeV to make an estimate of the decay widths. The 
values presented in Eq. (18) refer only to the last two columns in Table I that include the effect 
of the pion-exchange form factor. Thus, and considering both sets of errors, we cover the range 
0.5-2.8 MeV, which accounts for all values including errors quoted in the table. 
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FIG. 4. Dependence of the I 0 (m,D*o, m^o; Mx 2 , q^[mF 1 = ruF 2 = itid°]) integrand [see Eqs. (C5)-(C7)] on the 
pion loop momentum \l |. For the X 2 mass we have used 4013 MeV. Results with and without the inclusion 
of the pion form factor [Eq. (C7)] squared are presented in the right and left plots, respectively. In both cases, 
three different choices of the Gaussian regulator [Eq. (C6)] in the D*D*X 2 vertex have been considered. 

2. X 2 (4013) ->• DD*(D*D) 

Here, we will study the D + D*~, and D°D*° channels, which proceeds through the Feynman 
diagrams depicted in Fig. 5. This is also a d -wave decay so that both angular momentum and parity 
are conserved. The decay widths are expected to be comparable to those found for the A"2(4013) —>• 
DD decays, despite the phase space is considerably more reduced. This extra enhancement is caused 
by the extra multiplicity due to the spin of the final D*(D*) meson. 



FIG. 5. Feynman diagrams for the X 2 (4013) —> D + D* decay. Diagrams for the X 2 (4013) —> D°D*° 
transition are similar, with the appropriate changes of the exchanged pion charges. 

As commented before, we treat charm mesons non-relativistically and obtain the decay amplitude 
for the A2(4013) D + D*~ process as 

N*n 2 r 

T{ A,A*) d+d „_ = i—^e i j(\)e mn je* n (\*)lgc 2 P rn (m D ,+ ,m n o]Mx 2 ,q^) 

J 7T 

+2gX 2 I im (m D *o,m 7V --M X2 ,q»)} , (19) 

in the resonance rest frame. Here, q is the 4-momenta of the D + meson, and e n (A*) is the polarization 
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vector of the final D* meson with helicity A*, N* = 8M X2 m D* \Jmomf]* and ej mn is the 3- 
dimensional Levi-Civita antisymmetric tensor. Analogously, the X2(4013) —>■ D°D*° amplitude 


is, 


AT* q2 s 

T{\,\*) d o d *o = i^e ij (\)e mnj e* n (\, f ){2g^r m (rn D , + ,m n+] M X2 ,q^) 
J 7r 

+go 2 I im (m D * 0 , m n 0 ; M X2 

The two-body decay width in the rest-franre in this case reads [5]: 

dT<1 ~ jr^|T(A,A*)/ - — - a = D + D*-,D°D*° 


( 20 ) 


dQ(q) 


A,A, 


327t 2 M£ 2 ’ 


( 21 ) 


The sum over the D* and X 2 polarizations can be easily done, and we get 

^ ^ e^j(A)e rs (A) ^ ^ e (A*)e^(A *)tmnj£ips v T 2,5n5 rrn 3— 1,2,3 . 

A A* 

( 22 ) 

The above tensor structure should be contracted with I l7n (q) I rl (q). We see that the sum over 
polarizations of Eq. (22) is orthogonal to 5i m and 5 r i , which guarantees also here that only the UV 
finite Iq term of the three-point loop function is relevant. The contraction leads to a q 4 factor, 14 
and thus the width scales like \q | 5 , as expected for a d -wave decay. 

Results for the X 2 —» DD* decay widths are also compiled in Table I. We only show predictions 
for the ^(4013) —> D + D*~ (D°D*°) decays, because being the X 2 an even C-parity state, its decay 
modes into the charge conjugated final states have the same decay widths. In what respects to the 
effect of the form factors, the discussion runs in parallel to that in the Sect. IIIA1, though the 
effect of the pion-exchange form factor is significantly smaller here (a factor 4 or 5 at most). As 
expected, the widths are comparable to those found for the A 2 —>■ DD decays. Finally, we estimate 
the partial X 2 —> DD* (D*D) width (including both the charged and neutral channels as well as 
the charge-conjugated modes) to be 

T(A 2 -> DD*) + V{X 2 -> D*D) = (2.9 ± ^ j£8) MeV, (24) 

sys (A) 

where the errors have been estimated as in Eq. (18). The above result, together with that obtained 
previously for the DD channel, leads to a total X 2 width of the order of 2-8 MeV, assuming its 
existence. 


14 In the DD mode, studied in Sect. Ill A 1 a factor of 2<f 4 /3 is obtained instead. Thus, neglecting the D — D* mass 
difference, because the loop integrals are the same, we would find 


T t 


DD*(D*fl) 


| 2 ^fl Ton? 


(23) 


This extra factor 3/2 due to the spin-1 polarization vector produces an enhancement of the DD* decay mode with 
respect to the DD one, which partially compensates the smaller available phase space. 
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FFttI 

FFtt2 

FFttI 

FFtt2 

FFttI 

FFtt2 

FFttI 

FFtt2 


FFX 2 -G 

ffv 2 -g 

ffa 2 -l 

ffv 2 -l 

ffx 2 -g 

ffx 2 -g 

ffa 2 -l 

ffx 2 -l 

A [GeV] 

0.5 

0.5 

0.4 

0.4 

1 

1 

0.8 

0.8 

T(X 2 -+ D+D~) 

0 t 5+ 0 - 5 
U '°— 0.2 

0.6 

0.5 

0.6 

0 - 8+2 

01.0 

0.9 

0.9 

T(A 2 -+ D°D°) 

Q 4 + 0.5 
u '^- 0.2 

0.5 

0.4 

0.4 

0 - 6+2 

0.7 

0.6 

0.6 

T(X 2 -+ D+D*~) 

0 7 + 0.6 

u -'- 0.3 

0.9 (1.2) 

0.7 

0.8 (1.1) 

l-O+'i 

1.2 (1.6) 

1.0 

1.1 (1.5) 

T(X 2 -+ D°D*°) 

O 

l + 

0 0 

to 05 

0.7 (0.9) 

0.5 

0.6 (0.8) 

0-7+12 

0.8 (1.2) 

0.7 

0.8 (1.1) 


TABLE II. A'2(4013) —> DD 1 DD* decay widths (in MeV units) using different UV regularization schemes for 
the D*D*X 2 vertex and pion-exchange form factors. FFA^—G and FFX2—L stand for the results obtained 
with Gaussian (Eq. (5)) and Lorentzian (Eq. (26)) regularized local interactions, respectively. On the other 
hand, the widths in the columns FF 71 T and FF7r2 were obtained inserting the pion form factor of Eq. (C7) and 
F(l 2 , A w ) = e~ l in each of the two D*D^ 7r vertices, respectively. In the latter case, we take = 1.2 
GeV, as determined in the QCDSR calculation of Ref. [49] for the D*Dn coupling. In the DD* decay mode, 
we also show results (in brackets) obtained when a larger cutoff, = 1.85 GeV, is used in the D*D*tt 
vertex, as obtained in the QCDSR study carried out in Ref. [50] for this coupling. The results presented in 
Table I correspond to the [FFA2—G & FF7rl] columns and for the rest of choices we only provide central 
values. 


3. Charm decays: further analysis of uncertainties 

The uncertainties on the results compiled in Table I account both for reasonable estimates of 
HQSS corrections, as well as for the statistical errors on the inputs used to fix the LEC’s that 
determine the properties (mass and D*D* couplings) of the X 2 resonance. Moreover we are using 
an EFT to describe these decays, which means that there is an intrinsic uncertainty that can be 
determined systematically. We indicated the size of this error in Eqs. (18) and (24). For obtaining 
the EFT uncertainty we combined the predictions obtained for two different UV Gaussian cutoffs 
(A = 0.5 and 1 GeV) in the D*D*X 2 vertex, and used the spread of values to quantify this error. 
The idea is that the residual dependence of the results on the cutoff should provide an insight into 
the size of sub-leading corrections. 

Now, we try to further test the robustness of the systematic errors quoted in Eqs. (18) and (24). 
To that end, we have examined: 

• The dependence of our results on the functional form of the UV regulators, both in the 
D*D*X 2 and D*D^n vertices. 

— We have a contact theory with a Gaussian regulator and a cutoff A between 0.5 and 1 
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GeV. This theory, though very simple, is the LO of an EFT expansion describing the 
low energy dynamics of heavy hadrons (see Ref. [36] for details). Within the EFT we 
can expand observable quantities as a power series of the type 


a=Y, a -. 

V 



(25) 


where p is the momenta of the hadrons and A m the scale at which hadrons stop behaving 
as point-like particles (about the p mass). A LO calculation only keeps the first term in 
the series above. Hence it should have a relative error of order (p/A m)- We stress that 
this is expected to be so regardless of the exact regulator employed (gaussian, lorentzian, 
etc.), provided that the cutoff is at least of the order of Am- The reason for that is that 
the calculations we show are renormalizable: once the counter-terms are fixed 15 , they 
only contain negative powers of the cutoff A when we expand on the large cutoff limit. 
Hence, the uncertainty in the calculation is of order (p/ A). By taking A of the order 
of A Mi observables are guaranteed to have a cutoff uncertainty of the order (p/Am), 
equivalent to the EFT uncertainty. 

There are several methods for making a more precise estimation of the EFT error: the 
one we use in Table I is to vary the cutoff around values of the order of Am (hence 
the choice of the 0.5 — 1 GeV cut-off window). EFT predictions for two different cutoffs 
differ by a factor of (p/A m) and that is why we chose this particular way of assessing 
the errors. 

Alternatively, one could use different regulators or form factors to assess the size of this 
error. This idea also gives a cross-check of the previous error estimates based on varying 
the cut-off. We have employed a different regulator to check that our former estimation 
of the EFT errors is correct and to show that the particular regulator employed is not 
important. We closely follow here the discussion of Sect. VII of Ref. [37] and consider 
a Lorentz form for the regulator 


(p\V\p') = C IX 


r a 2 1 

\ A2 1 

A 2 + p 2 

A 2 + p' 2 


(26) 


with two values of the cutoff, namely A = 0.4 and 0.8 GeV, which were obtained by 
multiplying the Gaussian cutoffs by a factor of \J2/t: [37]. The resulting 16 widths are 
presented in Table II and turn out to be rather insensitive to the form of the regulator 
(this is to be understood within the limits of the expected EFT uncertainty). 


15 See for instance the discussion of Eqs. (23) and (24) of Ref. [37]. 

16 With the Lorentzian regularized local potential, we re-obtain the counter-terms Cox and C\x from the X(3872) 
inputs, which are then used to find the mass of the X 2 resonance and its couplings to the D*°D*° and D*~D* + 
meson pairs. These physical quantities hardly change, because the X 2 is a very loosely bound state and its dynamics 
is very little sensitive to the details of the D* D* interaction at short distances. 
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— Next we have studied the dependence of the widths on the pion form factor that accounts 
for the off-shellness of the pion in the mechanisms depicted in Figs. 3 and 5. To that 
end we used the results from the QCD sum rule (QCDSR) calculations performed in 
Refs. [49, 50]. The first of these two works considers the D*Dtt vertex: it was found the 
form factor is harder if the off-shell meson is heavy, implying that the size of the vertex 
depends on the exchanged meson. This means that a heavy meson will see the vertex as 
point-like, whereas the pion will see its extension. The authors of Ref. [49] find an on- 
shell value g = 0.48 ± 0.05 (note the different definition used in this reference), around 
1-2 sigmas below the value of 0.57 used in this work. In addition, they adjust their 
results for off-shell pions to a Gaussian form e ;2//A % with U l the pion four momentum, 
and find = 1.2 GeV. This form-factor 17 , in the region of interest for this work, 
turns out to be in a good agreement with that used to obtain the results presented in 
Table I. This can be seen in the new results showed in Table II and obtained with this 
new pion off-shell form-factor (FF7 t2 ). Deviations from our previous estimates of the 
Af2(4013) —> DD,DD* decay widths are both much smaller than the (HQSS & exp) 
uncertainties quoted in Table I, and well comprised within the systematic uncertainty 
generated when the D*D*X 2 cutoff varies in the 0.5-1 GeV window. 

In the X 2 —>■ DD* decay, there also appears the D*D*n coupling in one of the vertices, 
see Fig. 5. The off-shell behavior of this vertex might differ from that of the D*Dir 
one. This coupling was studied using QCDSR in [50] where, translating the definition 
used therein to the one used here, it was found an on-shell value g = 0.56 ± 0.07 in 
good agreement with the HQSS expectations. The off-shell behavior was described by 
a Gaussian, as in the case of the D*Dir vertex, but with a significantly larger cutoff, 
A 7r = 1.85 GeV. This significant difference is somehow surprising from the HQSS point 
of view, and we should note that the QCDSR actual calculations in [50] were carried 
out for significantly larger values of l 2 > 4 GeV 2 than in the case of the D*Dir coupling 
analyzed in Ref. [49]. Nevertheless, we used this softer dependence for the D*D*ir 
vertex and re-computed the DD* widths. Results are shown in brackets in Table II. 
Changes are now larger, and in general are of the order of 50%, though they could be 
still accommodated within the HQSS and EFT uncertainties already considered in our 
original calculations. The large momenta of the external mesons, that can even exceed 


17 To use the form factor of Ref. [49] in the computation of the widths, we have approximated the pion four momentum 
-*2 

squared by — l , which is sufficiently accurate because the dominant part of the integration comes from regions where 
the two virtual D* and D* mesons are almost on-shell. In this region, the energy of both heavy light vector mesons 
is approximately Mx 2 /2 which coincides with that of the external heavy mesons, and hence 1° is much smaller than 
1*1 . 
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0.5 GeV, make it possible that the short distance details of the decay mechanisms could 
be relevant. This is the weakest point in our scheme. The reason is that the EFT 
uncertainty is expected to be 0.5 GeV/A m > 1/2, as the variations of the DD* widths 
in Table II seem to confirm, and the calculation is really on the limit of validity of this 
kind of description and should only be considered as a reasonable estimate. 


• The contribution of decay mechanisms driven by the exchange of shorter range mesons (heav¬ 
ier) than the pion. 

Since the momenta of the external charmed mesons can exceed 0.5 GeV, one might think 
that shorter range contributions such as the p and uj exchanges could be sizable, and even 
comparable to those of the diagrams depicted in Figs. 3 and 5 for the exchange of a pion. 
We will focus on the X 2 —>• DD decay mode, where the momenta of the external mesons is 
the largest and we will begin by studying the effects of the exchange of a p meson. If we 
use the phenomenological D*Dp Lagrangian given in Eq. (3e) of Ref. [51], we find that the 
amplitudes of this new contribution can be obtained from those driven by pion exchange, and 
given in Eqs. (13) and (14), by replacing m n 0 and rri n + by m p 0 and m p + in the loop integrals 
and 


El 

P 

J IT 


m D 2 

— 9d*d p 
m D » H 


(27) 


where we have neglected \q | 2 /m^ w terms, with q the c.m. three-momentum of the external 
D and D mesons. The coupling constant gn*Dp has been computed in various schemes [51 — 
55] (ordinary and light cone QCDSR, vector dominance model, SU(4), etc.) and it varies in 
the range [2.8 ± 0.1,4.3 ± 0.9] GeV -1 (see Table 5 of Ref. [55]). Taking gr>*Dp ~ 5 GeV^ 1 , 
in the highest part of the interval of calculated values, we would have gjj*r> p /(g 2 / f 2 ) ~ 2/3, 
while a direct calculation of the loop integrals shows that those calculated using the p mass 
instead of the pion mass are around a factor of two smaller. Altogether, this indicates that 
the absolute values of the p-exchange amplitudes are about a factor of three smaller than 
those driven by the pion exchange. If one uses go*Dp ~ 3 GeV^ 1 , now in the lowest part 
of the interval of values, the p contribution will be, at the level of amplitudes, around eight 
times smaller than those considered in the present work. In any case, these effects are smaller 
than the HQSS and EFT uncertainties quoted in Table I, and therefore it seems justified to 
neglect them. On the other hand, the w-exchange contributions are even smaller, around a 
factor of three, because this meson is neutral, and it cannot generate mechanisms where a 
light charged meson is exchanged. 
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To estimate the size of the p exchange contribution in the X 2 —>■ DD* decay mode, in addition 
to the D*Dp vertex, we have used the phenomenological D*D*p Lagrangian given in Eq. (If) 
of Ref. [51]. The amplitudes of this new mechanism can be obtained from those driven by 
pion exchange, and given in Eqs. (19) and (20), by replacing m^o and m n + by m p 0 and m p +, 
and P 2 &) —> {P-iix) + x\q\/l) in the computation of the loop integral Iq (see Eq. (C5) of 
Appendix Cl) and using the appropriate coupling constants 


f 2 

J 7 r 



9D*Dp 


9d*d*p 

mo* 


(28) 


where we have neglected as before the |g | 2 /m^ w terms. (At this point, the phase in Eq. (28) 
is arbitrary, because the phase convention used in [51] and in this work for the heavy meson 
fields is different.) The gr>* d* p coupling is set to 2.52 in [51] , while it is estimated to be 4.7T0.2 
in [53, 56] (QCDSR). On the other hand, there exist some partial cancellations in the loop 
integral in this case, and it turns out to be at least a factor of 5 smaller than that calculated 
for the pion driven contribution. Thus, and having in mind the previous discussion for the 
expected range of values of the go* Dp coupling, we conclude that the p— exchange contribution 
can be safely neglected in the X 2 —>• DD* decay mode. 


In view of the results discussed in this subsection, we should acknowledge that as a result of 
the contribution from highly virtual pions, which is certainly in the limit of applicability of the 
low-energy EFT, the hadronic decay widths bear a large systematic uncertainty. Nevertheless, we 
have given arguments to be reasonably convinced that the results quoted in Eqs. (18) and (24) 
provide sensible estimates for the A2(4013) —>• DD,DD* widths. 

In the next subsection, we will study these hadronic decays in the bottom sector. There, the 
considerations are parallel to those discussed here for the charm sector. 


B. Bottom decays 

Thanks to the heavy flavor symmetry, the results of the previous subsection can be trivially 
extended to the bottom sector. There, we have a robust prediction, even when HQSS uncertainties 
(20%) are taken into account, for the X {,2 resonance (see Fig. 2). Moreover, all sort of non-relativistic 
approximations adopted in the current scheme are now more suited, since the range of variation of 
the internal pion momentum in the loops is similar to that shown in Fig. 4 for the charm sector. 

On the other hand, as discussed in Sect. IIC, we do not expect any significant isospin breaking 
effects and the X^ resonance would be a pure 1 = 0 state, with equal coupling to its neutral and 
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charged components. For simplicity, we will also neglect the tiny difference between B° and B± 
masses, and we will use a common mass m B = (m B o + m B ±)/2 = 5279.42 MeV. Yet, for the pion 
mass that appears in the loop integral, we take the isospin averaged value m^ = {2m n ± + m n o)/3. 
Note that the relevant internal pion momentum is around 750 MeV. With all these approximations, 
we find 


r(X b2 BB) = 


r(X b2 -> BB*) = 


3g 4 (y 62 ) 2 'rn B rrv Rt 


fn 


M x h 


q | 5 {lo(m B *,m n ;M Xb2 ,q fJ '[m Fl = m F2 = m B 


9 g 4 (g Xb2 ) 2 m,BTn%* 
M x h 


|<?| 5 [Io(m B *,m 7T -,Mx b2 ,q fJ '[mF 1 = m Bl m F2 = m B * 


(29) 


(30) 


lOvr/^ ^ Ab2 

for any charge channel (B + B~, B°B°, B + B*~, B°B *°) or charge conjugation mode ( B*B ). Our 
results for these decay widths are presented in Table III. We notice in passing that following heavy 
flavor symmetry we use the same value of g = 0.570 ± 0.006 in the charm and bottom decays. 
It agrees very well with a recent lattice calculation with relativistic bottom quarks which gives 
g b = 0.569 ± 0.076 [ 7] (we have added the systematic and statistic errors in quadrature). Yet, 
lattice calculations with static heavy quarks tend to give smaller values, see Ref. [48] and references 
therein. For instance, the ALPHA Collaboration presented a precise computation with the result 
of goo = 0.492 ± 0.029 [48]. Thus we expect that the decay widths of Table III slightly overestimate 
the real ones. 

For the BB mode we find a pronounced dependence on the Gaussian cutoff A employed in the 
dynamical generation of the resonance. This is inherited from the strong dependence of the X b2 
mass on this UV cutoff, as discussed in Eq. (12), which affects the available phase space for the 
decay. With all these shortcomings, we expect a partial width in the 1-10 MeV range, when both 
charge modes are considered. 

In the BB* decay mode, the impact of the Gaussian regulator is even larger, because it turns 
out that for A = 1 GeV, the central value of the resonance mass Mx b2 = 10594126 MeV is located 
below the threshold m B + m B * ~ 10604 MeV. Thus, in that case, the decay will be forbidden. For 
A = 0.5 GeV, we estimate a width also in the 4-12 MeV range, when the four possible decay modes 
(B+B*~, B°B*°, B~B* + , B°B*°) are considered. 


IV. THE X 2 AND X b2 RADIATIVE DECAYS 

In this section, we will study the X 2 —> DD *7 and X b2 —>• BB *7 decays. The interaction of the 
photon with the s-wave heavy mesons contains two contributions that correspond to the magnetic 
couplings to the light and heavy quarks [57] (see also Appendix A3). Both terms are needed 
to understand the observed electromagnetic branching fractions of the D* + and D*° because a 
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without pion-exchange FF 

with pion-exchange FF 


A = 0.5 GeV A = 1 GeV 

A = 0.5 GeV 

A = 1 GeV 

T(V 62 -s- BB) [MeV] 

26.0Zt|'.3 8 -7 5 

4 4+ 01 

0 7 + 1 - 4 
u -‘- 0.6 

T(V b2 BB*) [MeV] 

7 1 + 3 - 4 

1 • 1 -3.7 

2 0 +0 - 9 
z - u -1.0 

- 


TABLE III. Xb 2 —> BB, BB* decay widths (here BB^ refers to any of the neutral or charged modes, but it 
is not the sum of both). The results are given for different treatments of the three-point loop function. The 
errors have been obtained using the Monte Carlo analysis explained in Table I, but now considering the X 52 
mass histograms displayed in Fig. 2 and the coupling given in Eq. (12). The decay width of the V b2 BB* 
mode is the same because of charge conjugation symmetry. 


cancellation between the two terms accounts for the very small width of the D* + relative to the 
D*° [5]. Actually, one finds [57, 58]: 


T(D*° -»• D° 7 ) 
T (D* + -> D + 7 ) 


3 mu *0 \ om c J ' 
3 ni£)*+ \ 6m c J ' 


Pi 

P2 


2 _ g 2 m K 

3 P 87 Tf 2 K 

9 2 m n 

87T/2 


g 2 m n 

Svr /r 


(31) 

(32) 


where is the photon energy, m c the charm quark mass and a ~ 1/137.036 the fine-structure 
constant. In the non-relativistic constituent quark model P = 1 /m q ~ 1/330 MeV 1 , where m q is 
the light constituent quark mass. Heavy meson chiral perturbation theory allows one to improve 
upon this approximation by including corrections from loops with light Goldstone bosons, which 
give OQm q ) corrections [57]. 

Using isospin symmetry to relate r(H*° —>• D° 7 r°) and r(U* + —> D 0 n + ), correcting by the 
slightly different available p -wave phase space in each of the two decays, and taking into account 
the experimental D*° and D* + widths and radiative branching fractions quoted by the PDG [5], 
we find: 


T(D*° -> D° 7 ) = (22.7 ± 2.6) keV, T(D* + -> D + j) = (1.33 ± 0.33) keV. (33) 

These values differ from those used in Ref. [58] because of the recent accurate BABAR measurement 
of the D* + decay width, mentioned in Sect. Ill, which is around 10% smaller than the previous 
CLEO one used in Ref. [58]. Fixing the charm quark mass to m c = 1.5 GeV, we fit the parameter 
P to the above experimental values and find /3 _1 = (293 ± 11) MeV. 

In what follows, we will study decays of the type X 2 —> PP* 7 , being P and P* pseudoscalar and 
vector heavy-light mesons, respectively. Let us define p±, p% and p 3 as the four vectors of the final 
photon, pseudoscalar and vector mesons, respectively. Besides, let us define the invariant masses 













23 


m ij = Pij = (Pi +Pj) 2 i which satisfy 

m \ 2 + m \ 3 + 7TT-23 = M_y 2 + m\ + m 2 + m§ = ilf} 2 + nip + m 2 P » . (34) 


Since, as we will see, the Feynman amplitudes depend only on the invariant masses m 2 2 and m 23 
of the final 7 P and PP* pairs, respectively, we can use the standard form for the Dalitz plot [5] 


dT = 


(2vr) 3 32M 3 


-|T| dml^dml 2 , 


(35) 




with \T\ the absolute value squared of the decay amplitude with the initial and final polarizations 
being averaged and summed, respectively. Thus, we readily obtain 


r = 


i i 


r -Mi 


[ 2 dm 23 

J(rr 


K m i2) max _o 

dm\ 2 \T\ , 


(^12 )mi n 


(36) 


( 27 t ) 3 32 Mx 2 J(mp+mp*) 2 

where for a given value of m 23 , the range of m\ 2 is determined by its values when pp is parallel or 
anti-parallel to p 7 [5] : 


(m? 2 ) max = (e; + E*p ) 2 - (p; - p* P ) 2 , 

(m? 2 ) min = (E* + E* P f - (p; +P* p ) 2 , (37) 

with Ep = (m 23 + m 2 P — m 2 P ,) /2m23 and F* = (M^- — m 23 )/2m23 the energies of the P meson and 
photon in the m 2 3 c.m. frame, respectively, and p* p the moduli of the corresponding 3-momenta. 

Because of parity conservation, this is a p- wave decay and hence the photon momentum appears 
always in the amplitudes. In the X 2 rest frame, it is given by |p 7 | = P 7 = M — m 2 3 /(2Mx 2 )- 


A. X 2 (4013) -> DD *7 

We will first consider the ^(4013) — > D°D*° 7 decay, which proceeds according to the Feynman 
diagrams depicted in Fig. 6 . This decay can take place directly through the radiative transition of 
the constituent D*° as shown in Fig. 6 (a), which is the tree level approximation. However, there are 
other mechanisms driven by the DD* FSI. After emitting the photon, the vector meson D*° transits 
into the D°, and it can interact with the other constituent in the X 2 as shown in Fig. 6 (b). There is 
a third (c) mechanism in which the photon is emitted from the D*° meson, and the virtual D*°D° 
rescatter into D°D*°. Finally, Fig. 6 (d) and (e) present another possibility, namely the decay can 
also occur through the charged D* + D*~ component of the X 2 resonance, and the virtual charged 
D + D*~ and D* + D~ pairs then rescatter into D°D*°. Because the X 2 has a well defined charge 
parity (+), the decay width into the charge conjugated mode D°D*° 7 is the same. The Feynman 
diagrams contributing to the D + D*~ 7 decay mode are similar with obvious replacements. 
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FIG. 6 . Feynman diagrams for the X 2 (4013) —> D°D *°7 decay. Diagrams for the D + D* "7 transition are 
similar. 


One could also construct other FSI mechanisms by replacing the D*Dj vertices in Fig. 6 with 
the D*D *7 ones. These contributions are small and will be discussed below. 


1. Tree Level Approximation 


The Feynman amplitude of the mechanism depicted in Fig. 6 (a) reads (as in the previous sections, 
we treat the charm mesons non-relativistically) 


i'T (A, A*, 


^7 )d°d* 


gQ 2 {mi2)VknaN 1 



e ijm ei n (\)e* n (\*)e*j(\ 7 )p™ 
2 m D *o (mi 2 — m D *o + ie) 


(38) 


with mi 2 the invariant mass of the final 7 D° pair. Besides, e*(A 7 ) is the polarization vector of the 
final photon with helicity A 7 , p 7 is its three momentum and 1V 7 = 8 Mx 2 m %* \J m d^d* accounts 
for the normalization of the heavy meson fields and the Af 2 -D*°l )* 0 coupling. Finally, 


g 0 2 (mn) = go 2 x e -(Pl 2 2-7 2 )/A 2 = g x 2 x e -m D , 0 {m D » 0 -m 12 )/X 


(39) 


with p‘y 2 ~ rri D *o(Mx 2 — rn D *o — m\ 2 ) and y 2 = m D .o(My 2 -2m D .o) <0. The Gaussian form factor 
is inherited from the EFT UV renormalization scheme. 

We have neglected the D*° width in the above propagator because, since it is quite small, its 
inclusion only leads to small numerical variations in the decay rate which are certainly smaller than 
uncertainties induced by the errors in the coupling g^ 2 and the mass of the Af 2 (4013) resonance. 
Similarly, the use of the non-relativistic D*° propagator instead of (rn \ 2 ~ 1 leads also 






















25 


to numerically negligible differences, as compared to the HQSS corrections. The sum over the D *°, 
7 and X 2 polarizations can be easily done, and we get 


\r\ 


16tt aMxn 


D°D*° 


% 2 {rni2 


1 + 3^ 


(mi 2 - m D *0 + ie) 


2 ^ 7 " 


(40) 


The amplitude for D + D* 7 decay is readily obtained from Eq. (38) making the obvious re¬ 


placements: g'y 2 g ^ 2 , j3i —>■ ($2 and (m D o,m D *o) —> (m D +, m D *+). Performing the phase space 

integration, we find at tree level (assuming the existence of the X 2 state) 


T(A 2 (4013) -> D°D *° 7 ) tree = 18t§ (16±|) keV, (41) 

r(X 2 (4013)^ J D + H*- 7 ) tr ee = 0.10^° (0.09+°;°®) keV, (42) 


where the values outside and inside the parentheses are obtained with A = 0.5 and 1 GeV, respec¬ 
tively. The errors account for the uncertainty, both in the mass of the X 2 state and in its couplings 
g ^ 2 . derived from the A(3872) input ( M x t 3872) and the ratio R X { 3872) of the decay amplitudes for 
the A(3872) —>• J/ipp and X(3872) —> J/ipu decays) and the HQSS corrections, as explained in the 
the caption of Table I. We have neglected additional uncertainties stemming from the error on (3 
(~ 3%), because it is totally uncorrelated to those discussed above, and it is much smaller than 
those affecting for instance the g^ 2 couplings. The neutral mode width is much larger than the 
charged one, thanks to the bigger magnetic D*D^y coupling and a larger available phase space. 

In analogy with the discussion of Eqs. (20) and (21) in Ref. [33] for the A(3872) — »• D°D°tt 0 decay, 
in the X 2 radiative processes the relative distance of the D*D* pair can be as large as allowed by 
the size of the AT 2 (4013) resonance, since the final state is produced by the one body decay of the D* 
meson instead of by a strong two body transition. Thus, the radiative DD *7 decays might provide 
details on the long-distance part of the resonance wave function. For instance, the dT/d\pQ*o\ 
[dr/d\p£>*-\] distribution is related to the X 2 (4013) wave-function 'S’ipjj* 0 ) [^(pb*-)] [3 ]• This is 
in sharp contrast to the DD and DD* decay modes studied in the previous sections, which turned 
out to be strongly sensitive to short distance dynamics of the resonance, as revealed by the notorious 
dependence on the UV form factors. 

However, all these considerations are affected by the DD* FSI effects to be discussed next. 


2. DD* FSI Effects 

To account for the FSI effects, we include in the analysis the DD* —> DD* and D*D —» DD* T- 
matrices, which are obtained by solving the LSE [Eq. (A10)] in coupled channels with the VooOziM 
potential given in Eq. (A 8 ). Some isospin symmetry breaking effects are taken into account because 
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the physical masses of the neutral (DD*) and charged ( D + D* ) channels are used in Eq. (A 10 ). 
The X 2 -A D°D*° 7 decay amplitude for the mechanisms depicted in Fig. 6(b) and (c) reads 

- zT(A, A*, A 7 )^+ c) = U + e ij - m e^(A)€* n (A*)e; i (A 

x 4 m D m D * Too^oo(m 2 3) J{m D *o,m D *o,m D o,p y ) , ( 43 ) 

where m 2 3 is the invariant mass of the final D°D*° pair, 

200^00(^23) = T D ojj* 0 ^ ) , DO fi*o(m 23 ) + T d *o d°^d° d*° ( m 23 ), ( 44 ) 

and the three-point loop integral function J (Mi, M 2 , M3,p 7 ) is given in the Appendix C 2 . The 
integral is convergent, however, for consistency it is evaluated using the same UV regularization 
scheme as that employed in the £>(*).£)(*) EFT. In sharp contrast with the hadronic decays studied 
above, the momenta involved in these integrals are rather low. 

On the other hand, we see that in the (b)+(c) contribution there appears the combination 
T D OD*o^,D 0 D* 0 { m 23) + T D *o£)0_tE,o£)*o{m23)- In the isospin limit, when the mass differences between 
the neutral ( D°D*°) and charged ( D + D*~) channels are neglected, we will find Too^-oo = (^c=-i + 
T^l 1 _ 1 )/ 2 . From Eqs. (A8) and (A 10 ), we find the C -parity odd isospin amplitudes, 18 

[Tcwi] {m 23) = Cjz + G D 5*(m 2 3), I = 0,1 , ( 45 ) 

with G D £), — G D of)*o — G d +£)*~ a common loop function. Note that the kernel of this LSE is fixed 
by the isoscalar ( Cqz) and isovector (C\z) C(charge conjugation) = — 1 terms of V D (*)q(*). This is 
a trivial consequence of the conservation of this symmetry, taking into account that the Al 2 and the 
photon are even and odd C-parity states, respectively. 

The (d) and (e) FSI contributions of Fig. 6 are similar, with obvious replacements. We find 

- *T(A, A*,A 7 )^+; } = VWV 7 (fo + e Jim e^(A)e-(A*)e«(A>™ 

x 14 m D m D *T + — >0 o(m 23 )} J(m D *+,m D *+,m D +,p^), ( 46 ) 

where 

T)->00(^23) = T D+D *-^, D o£)*o(m 23 ) + T D * +D -_ }fD ofi*o(m 23 ), ( 47 ) 

and in the isospin limit, we would have T )_5.00 = i T c=°-i ~ T c^-i)/ 2 - 


18 Here and for simplicity we do not explicitly write the on-shell UV Gaussian form factors [see Eqs.(AlO) and (A13)]. 
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For the V 2 (4013) —>• D + D* 7 decay, the FSI contribution is 


- iT(A, A*, = ^4^N^i jm ^ n (X)e* n (K)e; i (X 1 )p. 


x Am D mn D >« gf 2 ( /3 2 + 


3 m ( 


T +—>+_ (m 2 3 )J(rn D *+, m D *+ ,m D +, p 7 ) 


+ 9(f 2 ( Pi + 


3 m, 


too- 


--(m 23 )J(m D *o,m D *o,m D o,p 1 ) 


(48) 


with 

Tj, - >+ -(m 23 ) = T D + D *-_t D + D *-(m, 23 ) + T D *+ D -^ D + D *~ (771,23), ( 49 ) 

Too^ + -(m 2 3) = T D oi)*o^ D + D *-(m23) ± ^D* 0 zj°->Z)+D*-( m 23 ) = T4 ->00(^123)- ( 50 ) 


and 74_ h _( m 23 ) = 700 ^ 00 (^ 23 ) in the isospin limit. 

Taking into account isospin corrections, induced by the meson mass differences, all the needed 

T-matrices (Too-»oo>^4_ >00 and T 4 _ >4 _) can be calculated by solving the coupled channel LSE, 

Eq. (A10), with the VdMdM potentials of Eq. (A 8 ), as mentioned above. Thanks to the conservation 
of C-parity, the FSI corrections will depend only on Cqz and C\z- Indeed, one finds 


TOO-OO ->00 

. ^bo-H— - +4 — 




Cgz+Ciz Cgz~C\z 
2 2 

Cqz—Ciz Cqz+Ciz 
2 2 


-1 


G d o d *° 0 
0 Gd+d*~ , 


^a(E) , 

(51) 


with .Fa (E) = Diag (/^ eu (77), /^ h (77)), where the Gaussian form factors are defined after Eq. (A13). 

The 24,(10610) observed in Ref. [59] carries electric charge, and its neutral partner was also 
reported by the Belle Collaboration [60]. It lies within a few MeV of the BB* threshold and it is 
tempting to speculate about it as a hadronic molecule. Belle also reported the discovery of a second 
exotic electrically charged bottomonium state [59], 24,(10650) in the vicinity of the 77*77* threshold. 
For both the 24,(10610) and 24,(10650) states, J F = 1 + are favored from angular analyses. 

Within our scheme, we assume that the 24,(10610) resonance is an isovector (7777* + 77*77) /\/2 
s-wave bound state with J pc = l -1 [15]. Note, HQSS predicts the interaction of the 77*77* system 

with 1 = 1, J pc = 1 + quantum numbers to be identical to that of the BB* pair in the 24,(10610) 
sector. Thus, HQSS naturally explains [31] the approximate degeneracy of the 24,(10610) and 
24,(10650). Taking for the 24,(10610) binding energy (2.0±2.0) MeV [62], we could fix a third linear 
combination of the LECs that appear in the LO Lagrangian of Eq. (A 6 ) 


C\ Z = C ia - C 1B = -0.75i°;[° (-0.30i°;° 2 ) fm 2 


(52) 


for A = 0.5(1.0) GeV. Errors have been obtained with a procedure similar to that described in the 
discussion of Eq. ( 8 ) and used in the case of the V(3872). Assuming heavy quark flavor symmetry, 
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the above value of C\z could be also used in the charm sector, subject to corrections of the order of 
0 (Aqcd/^c) — 20% that we take into account. Therefore, we predict the existence of the isovector 
charmonium partners of the Zf,(10610) and the Z;,(10650), though as virtual states in the second 
Riemann sheet [15], which probably correspond to the recently observed charged charmonium-like 19 
Z c (3900) and Z c (4025) states [63-6 ]. These resonances lie close to the DD* and D*D* thresholds, 
respectively, while J p = 1 + quantum numbers are favored from some angular analyses. 

Due to the presence of the Z c (3900) close to threshold, one should expect the loop (FSI) mech¬ 
anisms depicted in Fig. 6 to be important since T^T 1 _ 1 must have a pole. However, the value of 
Cqz = (Cqa — Cob ) is still unknown. It is not determined by the inputs deduced from the A)(3872) 
and Z ft (10610) states used in the present analysis. Depending on the value of Coz , there can be 
an isoscalar J pc = l -1 DD* s- wave bound state or not. For instance, considering the case for 
A = 0.5 GeV and taking the central value for Coz ~ — 2.5 frn 2 one finds a bound state pole in the 
DD* system bound by around 10 MeV; if Coz ~ —1.5 fm 2 , there will be a DD* bound state almost 
at threshold; if the value of Coz is larger, there will be no bound state pole any more. Therefore, 
the information of Coz will be crucial in understanding the DD* system and other exotic systems 
related to it through heavy quark symmetries [15, 35]. Conversely, as we will see, the X 2 radiative 
decay width could be used to extract information on the fourth LEC, Coz , thanks to the FSI effects. 

To investigate the impact of the FSI, in Fig. 7 we show the dependence of the partial ^2(4013) —>• 
DD* 7 decay widths on Coz■ For comparison, the tree-level results are also shown in the same plots. 
As expected, for the decay into the D°D*° 7 channel, the FSI effects turn out to be important, 
and for some values of Coz , they dominate the decay width. The maximum effects of the FSI 
mechanisms approximately occur for values of Coz which give rise to an isoscalar l -1 DD* bound 
or virtual state close to threshold. One can see an apparent deviation from the tree-level results in 
this region. When Coz takes smaller values, the binding energy of the bound state increases and 
moves apart from threshold; when Coz takes larger values, the pole moves deeper into a non-physical 
RS and becomes a virtual state further from the threshold. In both situations, the FSI corrections 
turn out to be less important. On the other hand, the FSI corrections are always important in the 
D + D*~ 7 channel. This is because the tree level amplitude involves only the D*±D±^ magnetic 
coupling, while FSI brings in the neutral magnetic coupling, which is much larger than the former 
one. This is also the reason, besides phase space, why the tree level width is much larger in the 

19 The simple EFT scheme employed in this work does not allow for finite width resonances, only for virtual or bound 
states. The merit of that EFT is actually making a connection between the ^(10610/10650) and (3900/4025) 
resonances on the basis of heavy flavor symmetry, not so much predicting the exact location. In that simple theory, 
we treat the £(,(10610) as a stable bound state and that is the reason why Z c ’s are predicted to be virtual states. 
Had we used a more complete EFT that takes into account the finite width of the .£(,(10610) and includes next-to- 
leading corrections, then we would have predicted the Z c ’s as resonances and located them more accurately in the 
complex plane. 
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neutral mode than in the charged one, as commented above. 

Notice that in the above calculations, we did not include the contribution from the coupled- 
channel FSI D*D* —>• DD* which can come from replacing the D*D 7 vertices in Fig. 6 by the 
D*D*^ ones 20 . We have checked that this is a good approximation when the re-summation of 
the charmed meson scattering is switched off. This is partly because the loop integral defined in 
Eq. (C 8 ) takes a much larger absolute value for the considered diagrams than for those with the 
D*D* 7 vertices in most regions of the phase space. The only exception is when the photon has a very 
low energy so that the D*D* are almost on-shell. However, after integrating over the phase space, 
this region provides a small contribution to the partial decay width because of the suppression 
due to relative p- wave between the photon and the l " 1 DD* system. The re-summation would 
introduce a further complication because of the presence of the Z c (4025) which couples dominantly 
to the D*D* in the l " 1 channel. It contributes mainly to the region close to the pole position of 
the Z c ( 4025). Again, this corresponds to the region with the very low-energy photon and thus it is 
suppressed due to phase space. Therefore, we expect that the neglected contribution discussed here 
has little impact, and its numerical effects should be safely covered by the sizable HQSS uncertainty 
exhibited in the results. 

To better understand the dependence of the ^2(4013) —>• DD* 7 decay widths on Cqz in the 
bottom panel of Fig. 8 we show the pole positions of the odd C-parity DD* —>• DD* T-matrix as 
functions of this LEC and for an UV cutoff of 0.5 GeV. There are two coupled channels, the neutral 
one which has the lowest threshold, and the charged channel. As a consequence, there are three 
relevant RS’s (among four). The first of them [(++)] is the physical one, while the two non-physical 
RS’s are reached by changing the sign of the imaginary part of the momentum inside of the loop 

functions G in Eq. (51) of the first channel [(—|-)] or the momenta of both channels [(-)]. Each 

of these non-physical RS’s are continuously connected with the physical one on the real axis above 
the relevant thresholds. The solid red and dashed blue curves stand for poles in the (+- 1 -) and (—h) 
RS’s, respectively. For sufficiently large negative values of Cqz, there is a bound state (a pole below 
threshold and located in the physical RS), which becomes less bound when \Cqz\ decreases. For a 
value of this LEC around —1.5 frn 2 , marked in Fig. 8 with a vertical black line, this state reaches 
the threshold. When Cqz is further increased, the pole jumps into the (—h) RS, becoming thus a 
virtual state (a pole below threshold and located in an non-physical RS) and it moves away from 
the threshold. In the medium panel we show the central values of the ^2(4013) radiative widths for 
both decay modes as a function of Coz, which were already presented in Fig. 7. As can be seen, it 

20 The eiectric part of the D* D* 7 vertex does not contribute to the FSI X 2 decay width amplitude when the quantum 
numbers of the final D* D* pair are T* , with the two heavy mesons in relative s— wave. Thus, we are left with the 
contribution from the magnetic coupling, as in the D*Dy case. 
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FIG. 7. Dependence of the X2(4013) —> D°D *°7 and the X2(4013) —> D + D *~7 partial decay widths on 
the low-energy constant Coz- The UV cutoff is set to A = 0.5 GeV (1 GeV) in the top (bottom) panels. 
The red error bands contain the DD* FSI effects, while the three horizontal blue lines stand for the tree 
level predictions of Eqs. (41) and (42). Besides the uncertainties on the mass and the couplings of the X 2 
resonance, the errors on C\z quoted in Eq. (52), together with the expected 20% heavy quark flavor symmetry 
corrections when it is used in the charm sector, are also accounted in the 68 % CL bands displayed in the 
panels. The red dash-dotted (full calculation) and solid blue (tree level) lines stand for the results obtained 
with the central values of the parameters. 

is around this critical value Coz = —1.5 fm 2 , when the FSI effects are larger for both decays, due 
to the vicinity of the pole to the threshold. This happens regardless of whether it is a bound or a 
virtual state, since the presence of the pole in both situations greatly enhances the odd C-parity 
DD* —>• DD* T-matrix near threshold. This can be appreciated in the top panel of Fig. 8 , where 
the dependence of the Too-s-oo and Tq_ *4 _scattering lengths on Coz is shown. 

Thus we have understood why in the region of values of Coz around —1.5 fm 2 , FSI corrections 
strongly affect the D°D *°7 decay width: this channel has the lowest threshold and the bound or 
virtual state is located on or nearby it. For the charged decay mode, the width exhibits a maximum 
for values of Coz also in this region, followed by a clear minimum placed now in the vicinity of 

Coz ~ C\z- Notice that when Coz = C\z, 7bo->q_vanishes (see Eq. (51)) and therefore the 

contribution due to the neutral mesons, driven by the largest magnetic coupling (/?i), in the FSI 
loops disappears. The exact position of the minimum is modulated by the further interference 
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FIG. 8 . Dependence on Cqz of several physical quantities predicted in this work. In all cases an UV cutoff 

A = 0.5 GeV is employed in the Gaussian regulator. Top: Too->.oo and _ >+ _ scattering lengths oat 

(red solid curve) and ac (blue dashed and green dash-dotted curves), respectively. They are defined as 
ai = HiTi(E = Mu + M 2 i)/ 2 jT, with /z, the corresponding reduced mass and (Mu, M 2 i) = (mco,mc*o) 
and (rri£)±,m£,t±) for the neutral and charged channels, respectively. The scattering length ac is complex 
because the neutral threshold is lower than the charged one, and therefore it is open. Middle: Central values 
of the X 2 (4013) — > D°D*° 7 (solid curve) and X 2 (4013) —> D + D*~"f (dashed curve) partial decay widths 
including FSI effects. Bottom: Position of the poles of the odd C-parity DD* —>■ DD* T-matrix, with respect 
to the neutral (me 0 + m^.o) threshold. Poles found in the various RS’s are shown (see the text for details). 
The red solid curve shows the evolution of the bound state with Cqz , while the dashed and the dash-dotted 
curves show that of the virtual ones. The vertical black line marks the value of Cqz for which a DD* bound 
state is generated at the D°D*° threshold. 

between the tree level and the FSI charged loops, which are comparable. 

In the bottom panel, we also observe a virtual state pole, the position of which is rather insensi¬ 
tive 21 to Cqz- It is originated by the interaction in the 1 = 1 sector, C\z, and it should be related to 

21 The situation is more complicated, as can be seen in the plot. There is a narrow region of values of Coz around 
[—0.4, —0.2] fm 2 , where the virtual state moves quickly away from threshold, shortly reappearing again (orange dash- 
dotted line) in a position similar to the one that it had at the left of Coz = —0.4 fm 2 . More in detail, in a narrow 

region included in the above interval of values of Coz, two poles (virtual) in the (-) RS coexist. Among these two 

poles, the decay width should be influenced only by the one closest to threshold, since it overshadows the second 

one placed deeper in the real axis. On the other hand, for Coz = —0.75 fm 2 the (-) RS virtual state (magenta, B) 

coincides with the one located in the (—h) RS (blue). This is easy to understand since at this point Coz = Ciz = Cz 
and then according to Eq. (51) the off-diagonal interaction term vanishes. In this situation, neutral and charged 
channels decouple, the scattering length ac is real, and the determinant (1 — CzG D op»o)(l — CzG D +p*~) would 
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the Z c ( 3900) exotic charmonium-like state reported by the BESIII and Belle collaborations. Within 
our LO EFT scheme, we do not find a DD* bound state, but instead a pole located near threshold 
in a non-physical RS [15]. 


vanish when either of the two factors in brackets is zero. It turns out that the factor (1 — CzG D o^.o) vanishes 

at the same energy both in the (—|-) and (-) RS’s, since by construction the G d oq,o loop function is identical 

in both unphysical sheets. However, the charged factor (1 — CzG D +£,,~) does not lead to any further pole for 
Coz = Ciz = —0.75 fm 2 and A = 0.5 GeV. In the A = 1 GeV case, not shown in Fig. 8, it happens that for 

Coz = Ciz = —0.3 fm 2 , there exist identical poles in the (—1-) and (-), and (-|—) and (-) RS’s, respectively, 

whose origin can be traced to the above discussion having in mind that now both terms in the decomposition of 
the determinant lead to poles. 
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B. X b2 ->■ BB *7 

The expressions found in the charm sector can be readily used here, having in mind the following 
correspondence: D° -H- B ~, D + -H- B°, D° B + , and D~ -H- f?°. Since the heavy quark flavor 
symmetry ensures that g and /?i i2 are the same in the b and c systems (up to corrections of the 
order Aqcd /me), the expressions of Eqs. (31) and (32) can be used to predict the widths for the 
B* radiative decays [57], 


r (B*~ £T 7 ) = (p 1 - -M E 3 = (0.49 ± 0.05) keV, (53) 

3 ms* \ 3m b J ' 

T{B*° -> S° 7 ) = f/? 2 - ^—) E ^ = (0.23 ± 0.02) keV, (54) 

6 ms* \ orrib) 


where we have taken the value = 4.8 GeV for the bottom quark mass. 

As in the study of its hadronic decays, we assume the Af b2 be a pure 1 = 0 state, with equal 
coupling to its neutral and charged components, = g Xb2 = -Eg Xb2 . The isospin breaking effects 
for the B* mesons are expected to be small and the tiny difference between the B° and B± masses 
can safely be neglected as well. In this limit we find at tree level 


/'/ (A, A*, A-y) 


DBB*^ 


( A 


1 A e ijm ei n ( A)e^(A,)ey(A 7 )py 
3 mb) 2rriB* {mw — ms* + is) 


(55) 


where f3 a = /3i(/3 2 ) for the B B* + 'y(B°B* 0 ^) mode, 8My b2 y/msms*, and mi 2 is the 

invariant mass of the final 7 f? pair. These amplitudes lead to 


T(X b 2 -> B~B*+^) tree = 13±H eV > r ( X b 2 -> B°5*° 7 ) tree = 6^° eV, (56) 


where the values have been obtained with A = 0.5 GeV. We remind that for A = 1 GeV, the central 
value of the resonance mass Mx b2 is located below the threshold (ms + rriB*) ~ 10604 MeV and the 
decay is forbidden. The errors reflect the uncertainty in the inputs from the V(3872) and the HQSS 
breaking corrections, as outlined in the caption of Table III; they are quite large and are dominated 
by those quoted for g Xb2 in Eq. (12). The widths are of the order of a few eV, significantly smaller 
than those of the radiative decays of the B*° and B*~ meson because of the quite reduced phase 
space available (~ 20 MeV) for this p -wave decay. They are also orders of magnitude smaller than 
T(V 2 -> D 0 D* 0 7 ) as a result of T(D*° -»• T>° 7 ) > T(B* -»• 5 7 ). 
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Coz (fm 2 ) 



Coz (fm 2 ) 


FIG. 9. Dependence of the Xb 2 —> B~B * + 7 and X^i B°B *°7 partial decay widths on the low-energy 
constant Cqz- The red error bands contain the BB* FSI effects, while the three horizontal blue lines stand 
for the tree level predictions of Eq. (56). 


The amplitude for the FSI mechanisms is readily evaluated and we find 


■Xb2 


iT (A, A*, A 7 )| S j* 7 = \/47raiV 7 ^=-ey m e : ' n (A)e* n (A*)e 7 (A 7 )p^J(m s *, m B *, m B ,p 7 ) 


x 4 m B m B * < ( Pi - 


3 mi, 


T c=- 1 (^ 23 ) ± T*d_i(rn 2 3 ) 


J BB* 


r c=-i( m 23) T Jc = Li(m 23 ) 

2 

where the -|— (—b) combination stands for the B~B* + 'y (B°B*° 7 ) decay mode and 77123 is now 
the invariant mass of the BB* pair. The C-parity odd isospin amplitudes are obtained by solving 
Eq. (45) using the bottom sector loop function G BB *. 

FSI corrections turn out to be important, as can be appreciated in Fig. 9. This is because we are 
generating in the T (7 ^ 1 _ 1 amplitude a bound state [Z&(10610)], almost at threshold (binding energy 
(2.0T2.0) MeV [62]), that enhances the loop mechanisms, as we discussed in the charm sector. If we 
pay attention for instance to the charged B~B* + 7 mode, we could appreciate a distinctive feature: 
there appears a destructive interference pattern between the tree level and the FSI amplitudes. 
Thanks to our MC procedure where correlations are consistently propagated, we also observe a 
reduction of the size in the uncertainties. Besides the uncertainties on the mass and the couplings 
of the Xb2 resonance, the errors on C\z quoted in Eq. (52) are also accounted for in the 68% CL 
bands displayed in the panels. Actually, these latter uncertainties should have also an important 
impact on the total CL bands. This is because variations of C\z allow for situations where the 
pole is located precisely at threshold (zero binding energy) or bound by about 4 MeV. In the first 
case the FSI contribution should be larger than that obtained with the central value of C\z, which 
correspond to a binding energy of 2 MeV. These big 68% CL bands makes hard to disentangle any 


+ ( & ~ i) 
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further dependence on Cqz, which in this case turns out to be quite mild. 


V. CONCLUSIONS 


In this work we have studied the hadronic and radiative decays of a molecular p*p* state with 
quantum numbers J pc = 2 ++ in the charm (X 2 ) and bottom (X^) sectors using an EFT approach. 
We have considered the V(3872) resonance as a J pc = 1 ++ DD* hadronic molecule. The X 2 and 
the Xb 2 states will be HQSFS partners of the V(3872) with masses and couplings to the P*P* 
heavy meson pair determined by the properties of the V(3872) resonance. 

The hadronic d- wave X 2 —>• DD and X 2 —>• DD* two-body decays are driven via one pion 
exchange. We observed that as a result of the contribution from highly virtual pions, which is out 
of control in the low-energy EFT, these hadronic decay widths (hence the total width of the X 2 
as well) bear a large systematic uncertainty. Even though the momenta involved in these decays 
probably lie outside the range of applicability of EFT the calculations are still valuable as a way 
to find reasonable estimates of these partial decay widths, which we expect to almost saturate the 
X 2 decay width. To this end and in analogy to the Bonn potential, we have included a monopole 
pion-exchange form factor, with a cutoff around 1 GeV, in each of the D*Dir and D*D*ir vertices 
to suppress the contribution of large momenta. We finally estimate the partial widths of both 
processes to be of the order of a few MeV. The analysis runs in parallel in the bottom sector with 
the assumption that the bare contact terms in the Lagrangian are independent of the heavy flavor. 
In this sector, we also find widths of the order of a few MeV. 

We discussed the radiative X 2 —> DD* 7 and Xb 2 —> BB *7 decays as well. The widths are small, 
of the order of keV’s (eV’s) in the charm (bottom) sectors. Furthermore, they are affected by the 
DD* or BB* FSI mechanisms. FSI effects are large because they are enhanced by the presence 
of the isovector Z c (3900) and 2^(10610) resonances located near the D°D*° and BB* thresholds, 
respectively. In the charm sector, FSI corrections turn out to be also sensitive to the negative C- 
parity isoscalar DD* interaction (Cqz)- Thus, future precise measurements of these radiative decay 
widths might provide valuable information on this LEC, which cannot be in principle determined 
from the properties of the V(3872), Z;,(10610) and 2^(10650) resonances. Constraints on this latter 
LEC are important in order to understand the dynamics of the p(*)_p(*) system. 
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Appendix A: Heavy meson Lagrangians: s-wave interactions and pionic and electromagnetic 

decays 

We collect in this appendix the Lagrangians used in this work. We use the matrix field 
[//(Q)] to describe the combined isospin doublet of pseudoscalar heavy-mesons P^ = ( Qu , Qd) 
[P(P' > = (uQ,dQ)\ fields and their vector HQSS partners P„* [P*,^] (see for example [40]), 



The matrix field H c [H c ] annihilates P [P] and P* [P*] mesons with a definite velocity v. The field 
Ha^ [Ha®' 1 ] transforms as a (2,2) [(2,2)] under the heavy spin <8 SU(2)y isospin symmetry [40], 
this is to say: 

H a Q) -> S (P (Q) ^) , H^ a -> (UH^Y St (A2) 

Their hermitian conjugate fields are defined by: 

H (Q)a = 7 °P,( Q)t 7 0 , Hp = 7°i7( < 3)“t 7 0 ) ( A3 ) 

and transform as [40]: 

H^ a -A (iJH^y S\ Hp -> S . (A4) 
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The definition for also specifies our convention for charge conjugation, which is CP^C 1 = 
p(Q) a anc [ CPa^C~ L = —P^' >a , and thus it follows 

CH^C- 1 = cH^ aT c-\ CH^ a C~ l = cffpc' 1 (A5) 

with c the Dirac space charge conjugation matrix satisfying cy^c -1 = — 7 J. 


1. Quadruple-heavy-meson contact interaction 

At very low energies, the interaction between a heavy and anti-heavy meson can be accurately 
described just in terms of a contact-range potential. The LO Lagrangian respecting HQSS reads [ 68 ] 


CiH = Ca Tr 

H { C) a H ( pC fii Tr 


+ C t a Tr 

V 

Tr 


+ C B TV 


Tr 

#(Q)a^(Q) 7 M 75 

+ CbTV 

7m75 

Tr \H^ )c T c d Hf ] ' 


(A 6 ) 


with r the Pauli matrices in isospin space, and C\ B light flavor independent LECs, which are also 
assumed to be heavy flavor independent. Note that in our normalization the heavy or anti-heavy 
meson fields, or H^\ have dimensions of E 3 / 2 (see [69] for details). This is because we use 

a non-relativistic normalization for the heavy mesons, which differs from the traditional relativistic 
one by a factor y/Mu- For later use, the four LECs that appear above are rewritten into Cqa , Cob 
and Cia, C\b which stand for the LECs in the isospin 1 = 0 and 1 = 1 channels, respectively. The 
relations read 


Co<j> — C<j> + 3CJ, Cirj, — Cfj, — CJ, for <p — A, B . (A7) 

The LO Lagrangian determines the contact interaction potential V = —£/4, which is then used as 
kernel of the two body elastic LSE (see Eq. (4) and the related discussion). 

The LECs that appear in the J P( = 1 ++ and 2 ++ sectors [Eqs. (3) and (9)] turn out to be 
Cox = Cqa + Cob and C\x = Cia + C\ B - The contact interaction in the Z;,(10610) sector (/ = 1 , 
JPC = 1+-) is C 1Z = C 1A - C 1B . 
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On the other hand, the interaction in the D*°D°, D + D* ,D* + D } space reads : 22 

UdMSO) = A T x Bia,g(C 0 z,C 0 x,Ciz,Cix) x A 


( Coa + C\a —Cob — Cib Cqa ~ C\a Cib — Cob ^ 

— Cob — Cib Cqa + C\a Cib — Cob Coa — C\a 

Cqa — C\a Cib — Cob Coa + C\a —Cob — Cib 

\ Cib — Cqb Cqa — C\a —Cob — Cib Cqa + C\a 

with Coz = Cqa — Cob and the orthogonal matrix A given by: 

(l 1 1 l\ 


(AS) 


A = 


(A9) 


1 1 - 11-1 
z 11 - 1-1 
V 1 -! - 1 l J 

Equation (A 8 ) trivially follows from the fact that the £ 4 h interaction of Eq. (A 6 ) is diagonal in the 
isospin basis and the charge conjugation is well defined . 23 The interaction given in Eq. (A 8 ) can be 
used as the kernel of an UV finite LSE to obtain the T-matrix that we use to account for the FSI 
in the radiative decays studied in Section IV, 


) D (*)] _1 - C{E 


[T dWSW (E)]- 1 =^ 1 (E)-{[V d( 
with the two particle regularized matrix propagator defined as 


^(E), 


a 


G(E) — Diag (G D o£,* 0 , G D *o D o t Gb+d*~ i G£,*+d-) > 
d 3 q e - 2 ? 2 /A 2 




(A10) 

(All) 

(A12) 


13 K ' 1 ( 2 vr ) 3 E - q 2 /2mj -Mi- Mj + ie ’ 

where trivially G D 0^*0 = G D *o^o and Gb+d*- = Gd*+d-- I n addition, the on-shell UV Gaussian 


form factor matrix reads 


E A (E) = Diag (fr(E),fr(E),ft(E)jt(E)) 


(A13) 


with f^ l \E) = exp(— k 2 /A 2 ) and k 2 = 2 fj, a (E — M\ a — M 20 ), with a = (neu), (ch). 


2 . pi.*)pi*) 7r interactions 


The relevant term in the LO Lagrangian of the heavy meson chiral perturbation theory [ 10-43] 
that provides the D*Dn and D*D*tt couplings is 


C„ HH = -JL(Tr 


1 + Tr [^(0)^(Q) 7 , 75 1 ) (? d ^) b a + • • • , (A14) 


22 In the bottom sector, the corresponding basis is: {B B * + , B* B + , B°B*°, B*°B 0 }. 

23 For instance in the charm sector, the C-parity states are [DD*]i : 2 = DD ^ D (1 ■<-»+, 2 — ). In our convention, 

the C-parity of these states is independent of the isospin and it is equal to =pl. 
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with (f> a relativistic field that describes the pion, 21 g is the heavy flavor independent PP*ir coupling 
and /,j- = 92.2 MeV the pion decay constant. Note that in our normalization, the pion field has a 
dimension of energy, while the heavy meson or antimeson fields or have dimensions of 

E 3 / 2 , as we already mentioned. 


3 . HH 7 interactions 


The magnetic coupling of the photon to the s-wave heavy mesons is described by the La- 
grangian [57, 58] 


£,HH = fQ a b [Tr 


+ 


eQ' 


4m, 


<2 


Tr 


}{{Q) b ]j(Q) 
fj{Q) a a n u fj{Q) 


+ Tr 
+ Tr 


jj(Q)b jj{Q) 


± fJLV 

Ffiu + • • • 


(A15) 


where = d^A u — d u A^ 1 with A M the photon field (CA^(x)C” 1 = — A^{x)), Q = Diag(2/3, —1/3) 
is the light quark charge matrix, and Q' is the heavy quark electric charge (in units of the proton 
charge e = \ZAna). For charm (bottom) Q' = 2/3 ( Q' = —1/3). Besides, mg is the heavy quark 
mass and f3 is the parameter introduced in Ref. [57]. These two terms describe the magnetic 
coupling due to the light (preserves HQSS) and heavy quarks (suppressed by I/toq), respectively. 
Both terms are needed to understand the observed electromagnetic branching fractions of the D* + 
and D*° because a cancellation between the two terms accounts for the very small width of the D* + 
relative to the D*° [5]. 

In the non-relativistic constituent quark model f3 = 1 /m q ~ 1/330 MeV” 1 , where m q is the 
light constituent quark mass. Heavy meson chiral perturbation theory provides contributions from 
Goldstone boson loops, which give O(y^m^) corrections to the decay rates [57]. If these loop 
corrections are evaluated in an approximation where heavy hadron mass differences are neglected, 
the correction to the above formulas can be incorporated by making the following replacements [57] 


PQn PQn ~ 

8 KJk 

PQ22 PQ22 + 


g 2 m 77 
8tt/2 ’ 


(A16) 

(A17) 


with f K ~ 1.2f n . 
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We use a convention such that <f> = ^ creates a 7r from the vacuum or annihilates a 7r + , and the c f> z field 
creates or annihilates a 7r°. We adopt the usual convention C(r • = (r • </>) T . 
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D* 




FIG. 10. The Xi self-energy diagrams from the s-wave D*D* and d-wave DD , respectively. 

Appendix B: Validity of the perturbative treatment of the DD for the 


In this appendix, we will argue that the d -wave DD may be treated perturbatively in the 2 ++ 
system. Even though this was already discussed in Ref. [7], we have included here a new argument 
grounded on a different EFT to make a more compelling case on the smallness of this contribution 
to the X 2 mass. We will compare the power counting of the self-energy diagrams of the X 2 from 
the d -wave DD and the s-wave D*D* two-point loops, see Fig. 10. If the DD loop is suppressed 
in comparison with the D*D* one, it will validate the perturbative treatment of the DD. Because 
in our case the heavy mesons are non-relativistic, we can apply a velocity counting for the loops 
analogous to the power counting of the heavy meson loops in heavy quarkonium transitions [70, 71]. 


For the D*D* loop, the velocity counting of the self-energy reads as 

S D*D* ~ 9 S( y 2\1 = 9s v i (Bl) 

where gs denotes the value of the s-wave coupling of the X 2 to the D*D*, v denotes the velocity 
of the D* meson, v 5 is for the loop integral measure since the non-relativistic energy is counted as 
0(v 2 ), and l/(y 2 ) 2 accounts for the two non-relativistic propagators. 

Similarly, for the DD loop, denoting the velocity of the D meson by w, the velocity counting is 
given by 




5 4 

0 W W 2 'S 

~ 3077^2 =9dUJ > 


( u > 2 ) 5 


(B2) 


where gp is the d -wave coupling constant normalized to have the same dimension as gs, and the 
factor of re 4 in the denominator comes from the two d -wave vertices. 


Therefore, we obtain the ratio 


r D/s = 


^DD 

^ D*D * 


g%w^ 

9 s v ' 


(B3) 


The question is now how gs> compares with gs- We can estimate go by considering the one-pion 
exchange diagram considered in this work as illustrated in Fig. 11. Because the X 2 is very close to 
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FIG. 11. The contact term of the d -wave coupling of the X 2 to the DD may be estimated by the one-pion 
exchange diagram. 

the D*D* threshold, we should count each of the D*(D*) propagators as 1/v 2 . This is equivalent 
to affirm that the cut due to the D*D* in the triangle diagram in Fig. 11 is the same as that in the 
D*D* bubble diagram of the X 2 self-energy. Thus, we can count the D* in both diagrams in the 
same way. But the pion propagator should be counted differently. The reason is that because the 
X 2 couples to the DD in a d- wave, the momenta in the D*Dir vertices of the one-pion exchange 
diagram should become the momenta of the D and D, qd = m^w, and the pion momentum is 
of the same order as we discussed in Sect. IIIA1. This is to say that the pion propagator should 
be counted as 1 /w 2 rather than 1/v 2 . Thus, expressing the content of Fig. 11 in terms of power 
counting gives 


9DW 


9S 


(v 2 ) 2 w 2 




(B4) 


where g is the axial coupling constant in Eq. (A14), and A a = 4 Tif n is the hard scale for the chiral 
expansion. 

Numerically, for the case of the X 2 , we have vj — yf ( Mx 2 — 2 mo)/mo — 0.38, and v — 
y/(M X2 - 2m D *)/m D * 0.06 if we take 7 MeV as the binding energy (recall that we have the 

charged D*D* channel explicitly whose threshold is around 7 MeV above the neutral one). With 
these values, we use Eqs. (B3) and (B4), which leads to go ~ 0.6 gg, to obtain an estimate of the 
contribution of the d -wave DD to the X 2 self-energy relative to the s-wave D*D*, 


r D /s ~ 0.05. (B5) 

The above value suggests a high suppression of the d -wave DD in comparison with the s-wave 
D*D*. We notice that the power counting of Ref. [7] indicates that the size of the DD loop 
is N 4 LO (next-to-next-to-next-to-next-to-leading order), in line with the velocity power counting 
arguments. 
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Appendix C: Three-point loop functions 
1. Hadron decays 


In this section we address the tensor three-point loop function that appears in the hadron 
decay amplitudes studied in the Sect. III. It is composed by a pion, and two heavy meson ( p*p *) 
propagators. The integral reads ( q = —k, q° + k° = Mj 2 ) 

d A l r V 


I 1 ' 3 [M. rri] Mx 2 , q >l ) = i 


■ / 


(2vr) 4 [(^ + q) 2 — M 2 + ie] [(k — l) 2 — M 2 + ie] ( l 2 — m 2 + ie) 

_ tl _ (on 

4 M 2 J ( 27 r) 4 (i° + q° — uJh + is) (k° — /° — uJh + is) (l 2 — m 2 + ie) ’ 
where M is the mass of the heavy particles in the loop, m is the mass for the light intermediate 
particle, Mx 2 is the total c.m. energy, and q and k are the external 4-monrenta of the two particles 
in the final state of masses and mp 2 , respectively. In addition, 


q o = Mx 2 + m 2 Fl - m F2 


k o = M 2 x 2 + rn% - m 2 Fi 


(C2) 


2 Mx 2 2 Mx 2 

and Uh = M + (q + l) 2 /2M is the non-relativistic energy of the virtual heavy mesons. Using 
Cauchy’s theorem to integrate over the virtual pion energy 1°, we obtain : 25 




qv ) - — [ 

,q 1 ~ AM 2 J 


d 3 l 


ri j - 


Mx 2 — 2 uj h — 2w 


(C3) 


(27t) 3 2 Lo(k° - u> - u h )(q° - uj - u h )(M X2 - 2u h ) 

for \m 2 Fi — rripj < 2 mMx 2 < AMm to guarantee that the integral in Eq. (Cl) is real, and ui(l ) = 
y/m 2 + l 2 . The loop integral I i3 presents a logarithmic UV divergence. Indeed, P 3 admits a tensor 
decomposition 


I l3 {q) = h{q 2 )q\ 3 + h{q 2 )d 13 |g| 


(C4) 


The I\ term presents an UV divergence, but it does not contribute to the amplitude because it 
annihilates the traceless spin-2 polarization tensor. This means that only the Iq term is relevant. 
It can be computed as: 


I 0 {M,m;Mx 2 ,q 2 ) 

^ 1 [ + °° dll* f +1 (x) _ M Xa -2u h (l,x)-2u _ 

32M 2 7r 2 q 2 J 0 u ; J_! ^ (k° - uj - uj h (l,x))(q° - u - u h (l,x))(M X2 - 2u h (l,x))' 1 

with ujh(l, x ) = M+(l 2 +q 2 + 2\l\\q \x)/2M, and P 2 the Legendre’s polynomial of order 2 . This term 

is not UV divergent because in the limit l —» +00 all dependence on x, besides P‘ 2 (x), disappears and 


If ttifi = m,F 2 , k° = q° = Mx 2 /2 and the loop function now reads, 


r j (M,m-,Mx 2 ,q) ^ 


b! 


Pi 


rv 


8 M 2 J (27r) 3 u}(Mx 2 /2 — uj - u h )(Mx 2 /2 — Uh) ’ 


Mx 2 < 2 M . 
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the integration over x gives zero. The convergence of the integral is greatly enhanced because P 2 (x) is 
orthogonal to x as well. Moreover, the same type of arguments guarantees that Iq(M, m\ Mx 2 , Q ' 2 )~ 
const, in the q 2 —>■ 0 limit. Numerically, we use non-relativistic kinematics to compute \q | in the 
evaluation of Iq{M, m; Mx 2 , Q 2 ) in Eq. (C5), i.e., q 2 — 2 /if 1 f 2 (Mx 2 — n r iF 1 — rriF 2 )- However, to 
guarantee the appropriate d- wave phase space, we use relativistic kinematics to evaluate q L q 3 in 
Eq. (C4) and the \q | phase-space factor that appears in Eqs. (15) and (21). 

For consistency with the scheme adopted in Eq. (5), we include a Gaussian regulator in the 
p*p*X 2 vertex by multiplying the integrand in Eq. (C5) by the exponential factor, 

e —(?+f) 2 /A 2 e -(P»+^+2|nif|*)/A» 

g— 7 2 /A 2 = g— 7 2 /A 2 

with 0 > y 2 = M(Mx 2 — 2 M). We divide by the factor e ~^ 2 / a2 , because it was incorporated in the 
p* p*X 2 coupling. 

In addition, the exchanged pion is highly virtual, and one might include a vertex form factor of 
the form 

F(f 2 , A,) = K ~ 1 GeV, (C7) 

l 2 + A 2 

in each of the two it p(*)p(*) vertices. 


2. Radiative decays 


In the computation of the FSI effects on the radiative decays of the X 2 and Xy 2 resonances in 
the Sect. IV, the following three-point loop function appeared 


J (Ml , M 2 , M 3 , Pry) 


d A q 1_1_1_ 

(27t) 4 q 2 — M 2 + ie (p — q ) 2 — M| + is (q — p 7 ) 2 — M| + ie 


(C8) 


i 

8M1M2M3 


d 4 q 1 1 

(27t) 4 q° - Mi - q 2 /2Mi + ie My 2 — q° - M 2 - q 2 /2M 2 + ie 


1 

q° - Ery- Ms - ((? - p 7 ) 2 /2M 3 + ie 

j»i2/^23 f d 3 q _1_ 

2MiM 2 M 3 J (2vr) 3 (g 2 + c - ie) (g 2 


1 

2 /J -23 g • P 7 /M 3 + d - ie) ’ 


(C9) 


with = (My 2 ,0) in the rest frame of the V 2 and /j ^ 1 = (M ? 1 + M- 4 ). In addition, i > 12 = 
Mi + M 2 - M.y 2 , 6 23 = M 2 + M 3 + E 7 — Mj 2 , c = 2pi 2 bi2 and d = 2/x 23 6 23 + p 23 p 2 /M 3 . Since 
all the intermediate mesons in the present case are non-relativistic, the three point loop has been 


treated non-relativistically. This loop integral is convergent and its analytic expression can be found 
in Eq. (A2) of Ref. [71]. However, for consistency, despite the three-point loop function in Eq. (C 8 ) 
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being finite, it should be evaluated using the same UV renormalization scheme as that employed in 
the D^*’D^ EFT. This is accomplished by including in the integrand of Eq. (C 8 ) a Gaussian form 
factor, F\(q) defined as 

F A (q) = e -(9 2 -7 2 )/ A 2 e -(9c 2 m-e 8 heii)/ A2 . (CIO) 


Here y 2 = 2 /h 2 (M a - 2 - M 1 - M 2 ), g 0 2 nshcll = 2 ^ 23(^23 - M 2 - M 3 ), with m| 3 
— 2Mx 2 E 7 , and 

^2 _ Pi) 2 + ^3g 2 

Qcm “ M 2 + M 3 


(*> - P7) 2 = 


(Cll) 


Note that the first exponential factor accounts for the off-shellness in the X 2 D*° D*° coupling, as in 
Eq. (39), while the second one accounts for the virtuality of the incoming mesons in the DD* —» DD* 
and D*D —>• DD* T-matrices. Note that, after the inclusion of this factors, an analytical expression 
for the integral cannot be easily obtained, and it needs to be computed numerically. 
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